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Abstract 

This article gives a proof of Fejes Toth's Dodecahedral conjecture: the volume 
of a Voronoi polyhedron in a three-dimensional packing of balls of unit radius is 
at least the volume of a regular dodecahedron of unit inradius. 

1 Introduction 




Figure 1 : Model of Voronoi cells by Michael Knauss and Silvan Oesterle. 

A packing of congruent unit radius balls in three-dimensional Euclidean space is 
determined by and is identified with the set A of centers of the balls. A packing A 
determines a region called the Voronoi cell around each ball. The Voronoi cell Q(A, v) 
around a ball at v e A consists of points of space that are closer to v than to any other 
w e A. The Voronoi cell is a convex polyhedron containing v. Figure [T] shows the 
Voronoi cells of a finite packing. 

The Dodecahedral conjecture asserts that in any packing of congruent balls of Eu- 
clidean space every Voronoi cell has volume at least that of a regular dodecahedron of 
inradius 1 (Figure |2j. This bound is realized by a finite packing Ad,,d (of twelve balls 
and a thirteenth at the origin) obtained by placing a point of A^ at the center of each 
face of a regular dodecahedron (of inradius 2). 

The assertion can then be stated as the inequality 

vol(Q(A, v)) > vol(Q(Arf od , 0)) 
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for every v e A, and for every set of points A c R 3 whose pairwise distances are at 
least the diameter 2. The case of equality occurs exactly when Q(A, v) is congruent to 
a regular dodecahedron of inradius 1 . 

1.1 History 

L. Fejes Toth made the conjecture in 1943 Il34l . In that article, L. Fejes Toth sketches 
a proof based on an unproved hypothesis. This hypothesis is a quantitative version 
of the kissing number problem in three dimensions. This unproved hypothesis is now 
generally regarded as being nearly as difficult as the Dodecahedral conjecture itself. 

L. Fejes Toth returned to the Dodecahedral conjecture in a number of publications. 
It is a prominent part of his two books |36], 1 35 1. According to the strategy of 11361 . 
the Dodecahedral conjecture forms a step towards the solution of the sphere packing 
problem (discussed below). In [35] , he proved that the Dodecahedral conjecture holds 
for every Voronoi cell with at most twelve faces. This result is reviewed in Section[4] 
It is an ingredient in the proof presented here. 



Figure 2: The regular dodecahedron Voronoi cell and enclosed sphere 

A lower bound on the volume X of a Voronoi cell implies an upper bounds on the 
density 4-n/(3X) of packings of congruent balls in three dimensions. The Dodecahedral 
conjecture gives an upper bound on density of 0.755. Upper bounds on the density 
based on lower bounds on the volume of a Voronoi cell in the literature include Rogers' 
upper bound 0.7797 ||33j, and Muder's upper bounds 0.77836 (271 and 0.7731 EST 

In 1993, Hsiang published what he claimed to be proofs of the Kepler conjecture 
and the Dodecahedral conjecture H22H . However, the proof did not hold up to careful 
analysis. 'As of this writing, Kepler's conjecture as well as the dodecahedral conjecture 
are still unproven" p761]. See also, IT2I . 

An alternative approach to the Dodecahedral conjecture is described in Q]. Unfor- 
tunately, a counterexample has been found to both parts of the third conjecture of that 
article. The counterexample is described in the preprint [ 16] 

K. Bezdek conjectures that the surface area of any Voronoi cell in a packing of unit 
balls is at least that of a regular dodecahedron of inradius 1. This strengthened version 
of the Dodecahedral problem is still open Q. 
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1.2 The sphere packing problem 



The Kepler conjecture, also known as the sphere packing problem, asserts that no pack- 
ing of congruent balls in three dimensions has density greater than the density of the 
face-centered cubic packing. S. McLaughlin carried out the research for the proof of 
the Dodecahedral conjecture at the University of Michigan while S. Ferguson and T. 
Hales worked on the sphere packing problem. Both problems were solved in 1998. 

There is no strict logical connection between the two problems. The Dodecahedral 
conjecture does not follow from the Kepler conjecture and is not an intermediate step 
in the solution to the Kepler conjecture. (In Fejes Toth's strategy, it was an intermedi- 
ate step; however, that strategy was not followed in the solution of the sphere packing 
problem.) Nevertheless, the two solutions follow a similar outline and share a signif- 
icant number of methods. Both are based on long computer calculations. Computer 
code was freely exchanged between the two projects. 

This article is written in a way that it is not necessary to read or understand the 
solution of the sphere packing problem before reading this article. However, for the 
benefit of the reader, this article points out parallels with the sphere packing problem. 
It also cites various results from that proof. 

1.3 Differences 

Although the proof of the Dodecahedral conjecture runs parallel to the solution to the 
sphere packing problem, there are special difficulties that arise in the proof of the Do- 
decahedral conjecture. In no sense is it a corollary of the sphere packing problem. In 
the packing problem, there turn out to be many ways to reduce the infinite ball prob- 
lem to a problem about finite clusters of balls. This multiplicity of choices makes it 
possible to design many difficulties away. If one reduction is not satisfactory one can 
work with another. With the Dodecahedral problem, there is no such flexibility. The 
problem about finite clusters of balls is fixed from the outset. This gives the problem a 
degree of rigidity that is not present in the sphere packing problem. 

1.4 Ten years later 

Over ten years have elapsed from the completion of the research until publication. A 
few words of explanation are in order. The review and publication process for the Ke- 
pler conjecture extended from 1998 until 2006. Because of significant sharing between 
the Kepler conjecture and the Dodecahedral conjecture, a "wait and see" attitude devel- 
oped toward the Dodecahedral conjecture. Once the Kepler conjecture was published, 
the path became clear for the publication of Dodecahedral conjecture. 

The details of the proof in the current publication are essentially the same as those 
in the preprint posted on the arXiv ifTTl in 1998. No significant errors have emerged in 
the original 1998 preprint. However, the article has undergone various major rewrites 
since then, including a significantly expanded version in 2002 [16]. That version is 
written in such a way that it is entirely independent of the proof of the Kepler con- 
jecture. As an abridgment of a longer version, this article replaces some proofs with 
summaries and references. This article refers the reader when necessary to relevant 
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passages in the full version of the proof. The computer code has also been entirely 
rewritten. 

A formalization project, called Flyspeck, aims to provide a complete formalization 
of the proof of the Kepler conjecture |[T5l [TP . (A formal proof is one in which ev- 
ery logical inference of the proof has been independently checked by computer, all the 
way to the primitive axioms at the foundations of mathematics.) A parallel project, 
called Flyspeck Light, aims to do the same for the proof of the Dodecahedral con- 
jecture. These long-term projects will take many years to complete. Nevertheless, 
significant progress has already been made toward the formal verification of the com- 
puter code [29 , 30 1 . The revisions in this article incorporate the parts of Flyspeck Light 
that have already been completed. 

1.5 Truncation 

The distance from the center of the regular dodecahedron to a vertex is tdod = V3 tan(7r/5) ~ 
1.258. This parameter is used to truncate Voronoi cells; it makes volumes easier to es- 
timate. A similar truncation takes place in the solution to the packing problem with 
truncation parameter to = 1 .255. It is a happy coincidence that these two truncation pa- 
rameters are so close to one another. A great deal of duplicated effort might have been 
avoided if these two parameters were equal. However, the parameter tdod cannot be 
replaced with anything smaller, and although the parameter to could easily have been 
made larger, its value was already too deeply entrenched in published articles by the 
time work started on the Dodecahedral conjecture. 

As a first step towards unifying the proofs, many results can be stated in a form 
that holds for all t e [to, td„d]- To transfer a lemma from ifTSl to this article, a simple 
process is involved. The first step is generalization, replacing the constant to with a free 
parameter t € [to, tdod]- The second is specialization, 1 1-» tdod- 

The number f ( >, although rational, can be consistently treated as an independent real 
transcendental in the solution to the sphere packing problem; that is, none of the proofs 
involving to rely on its exact numerical value. The constant fo can always be replaced 
by a constant t in a suitably small interval about fo- However, the only way to know that 
this small interval is wide enough to contain tdod is to study the details of the proof. We 
have made a detailed study of relevant proofs in [ 18 1 to insure that they can be adapted 
to the present situation. 

As a matter of terminology, a proposition for the Dodecahedral conjecture is said 
to be a reparametrization of a proposition with the constant to, if it is obtained by 
mechanically replacing to with tdod, wherever that constant appears, and if the proof 
goes through verbatim with this minor change. When this occurs, there is nothing to 
further to be learned by repeating the proof. The cited proposition already contains all 
the needed detail. 

1.6 Terminology 

Various notation and terminology is shared between the solution of the sphere packing 
problem and this article. Vocabulary can be imported from the sphere packing problem 
three different ways. The simplest way to import a term is for the term to have precisely 
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the same meaning in both places. For example, the terms "orientation" and "packing" 
have the same meaning in both places. 

The second way to import a term is by making a reparametrization of a term that 
depends on the parameter f(j. For example, the definition of standard component in this 
paper is the reparametrization of standard component in [18]. 

The third way in which terms have been imported into this proof from the proof 
of the Kepler conjecture has been by structural analogy. For example, the term tame 
graph is a technical notion that arises in the solution of the sphere packing problems. At 
the analogous point in this proof, a somewhat related collection of graphs appears. To 
emphasize the analogous role they play in the proof of the Dodecahedral conjecture, 
they are called tame Voronoi graphs. Similarly, there is a strong analogy between 
Voronoi weight assignments in this article, and admissible weight assignments in ifTSl . 

There are a some terms in that have been renamed in this version for greater preci- 
sion. Special has become unstable, standard region has become standard component, 
dihedral angle has become azimuth angle, and non-external has become internal. Pla- 
nar maps have been replaced with hypermaps. A few terms (such as distinguished) 
have been slightly redefined, when doing so is harmless. 

2 Outline 

This section gives a precise statement of the main theorem and describes the broad 
outline of the proof. 

2.1 Formulation 

This article proves the Dodecahedral conjecture in a stronger version than that stated 
in the abstract. A truncation of the Voronoi cell already has volume at least as great as 
that of the regular dodecahedron. This subsection describes the truncation and states 
the stronger version of the main theorem in a precise form. 

Let A be a packing and let vo e A. Let B(vq, r) be a closed ball of radius r centered 
at vo. Let A(vo,r) = A n B(vo,r). Let S(A, vo) be the set of all S = {vo, Vi, V2, V3} c 
A(vq, ltdod) consisting of four distinct points such that |vj - V;| < ltdod for all i, j and 
such that the circumradius of each triangle {v,, V;, v^} c S is at most V2. Write conv(5') 
for the convex hull of S e S(A, Vq). 

Define the following truncation Q. trunc {A, vo) of the Voronoi cell Q(A, vo): 

{x e Q.(A, v Q ) : (x e B(v , t dod )) V (3S e S(A, v ). x e conv(S)) 

That is, the the Voronoi cell is truncated by intersecting it with a ball of radius tdod, 
except inside regions protected by the sets S e S{A, vo). Note that the special packing 
A dod satisfies O fr „„ c (A^, 0) = Q.(A dod , 0). The Dodecahedral conjecture takes the 
following strengthened form. 

Theorem 2.1. For every packing A and every v s A, 

vol(O franc (A, v )) > vol(Q,(A dod , 0)). 
Equality holds exactly when £l(A, vq) is congruent to Q.(A dod , 0). 
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2.2 Proof outline 



The Lebesgue measure is translation invariant. Thus it does no harm to assume that the 
center point Vo of the Voronoi cell lies at the origin: vo = € A. The assumption that 
vo = e A remains in force for the rest of this article. 

The next reduction is to replace the set A with A(0, ltdod). This is accomplished by 
the following lemma, which shows that the volumes in Theorem 2.1 are insensitive to 



points of A outside A(0, ltdod)- The proof appears in Lemma 5.3.6 
Lemma 2.2.1. Let A be any packing (with e A). Then 

n trunc (A,0) = £W(A(0,2f dod ),0). 

The condition A = A(0, ltdod) is a standing assumption for the rest of the article. 
Let A* = A \ {0} and let n = #(A*) be the cardinality. The proof of Theorem [ZT] splits 
into two main cases: n < 12 and n > 13. In fact, L. Fejes Toth settles the case n < 12 
in his book [35 1. Section[4]sketches Fejes Toth's proof. 

Completely different methods treat the case when n > 13. This part of the proof is 
considerably more difficult than the case treated by Fejes Toth. Here is a sketch of the 
proof of the case n > 13. This rough sketch will be expanded in greater detail later in 
the article. 

Let A be a packing satisfying the standing assumptions that e A and A = 
A(0, ltdod)- The Dodecahedral conjecture seeks a minimum to the objective function 

vol(O frmc (A,0)). 

This is a nonlinear optimization problem in finitely many variables. The target value 
for the minimization is u>dod = v °l(^(A/od, 0)). When n > 13, this article prove^j] 

vol(Q, r „„ c (A, 0)) > a> dod . 

Some combinatorial information about each packing A is encoded as a graph. The 
vertex set of the graph is A*. (Because of this graph, elements of A* are generally 
called vertices.) The edge set is 

E = {{v, w] : v, we A*, \v - w\ < ltdod}- 

This is a planar graph. The Dodecahedral conjecture reduces to the case when this 
graph is connected. In fact, when the graph is not connected, this article constructs 
another packing A' whose graph is connected, with the same cardinality as A, and such 
that 

vol(£W(A,0)) = vol(£W.(A',0)). 

Similarly, the conjecture reduces to the case where the graph is biconnected. This 
again involves constructing an auxiliary packing A" of the same cardinality and whose 
truncated Voronoi cell has the same volume. Now assume that the graph of A is bicon- 
nected. 



1 An examination of the proof shows the right-hand side can be improved to oJdod + 10 10 . In fact, in terms 
of notation to be described later in the article, the proof only relies on the bound /j(A, Uf) > for triangles 



F, but in fact every configuration has at least one triangle F with /i(A, Uf) > 10 



-10 
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Figure 3: Sphere packing A^d, Voronoi cell £l(Ad d), and planar graph G(A^) 

This article supposes the existence of a counterexample A to the Dodecahedral 
conjecture and makes a detailed study of the properties of its graph. It defines a class 
of graphs (called tame Voronoi) and proves that the graph of every counterexample is 
tame Voronoi. 

Tame Voronoi graphs can be described in purely combinatorial terms, without ref- 
erence to packings, Voronoi cells, and volumes. All tame Voronoi graphs can be clas- 
sified up to isomorphism. This classification is one of the main steps of the proof. 
There are only finitely many possibilities. Thus, the graph of any counterexample to 
the Dodecahedral conjecture must be one of these finitely many cases. 

Each tame Voronoi graph can be encoded as a hypermap H. (A hypermap can be 
defined as a finite set, together with two permutations on that set. The elements of the 
given finite set are called darts.) If H is a hypermap, let V be the finite dimensional 
vector space of real-valued functions on its set of darts. The pair (H, <t>) is called a hy- 
permap system if O is a set of boolean valued functions (f> : V ( — > {true, false) for some 
I. A hypermap system (H, <£) is feasible if there is some x = (x\, , . . , X{) e V e such that 
<p(x) holds for all (5 6$. The computer code repository specifies a hypermap system, 
called the Voronoi hypermap system, for each case H arising in the classification of 
tame Voronoi graphs. 

Another major step of the proof is the proof that every Voronoi hypermap system is 



infeasible (Theorem 8. 1 1. The proof of this theorem is a case-by-case analysis based on 
the explicit enumeration of tame Voronoi graphs, up to isomorphism. The feasibility 
problem for each Voronoi hypermap system is converted to a system of linear programs. 
The infeasibility of the Voronoi hypermap system follows from the infeasibility of the 
corresponding linear program. 

If there exists a counterexample to the Dodecahedral conjecture, there is an associ- 
ated Voronoi hypermap system (H, Q>). By the preceding result, this hypermap system 
is infeasible. On the other hand, the counterexample can be used to construct a feasible 



solution to the system (Theorem 8.2 1. This contradiction shows that a counterexample 



cannot exist. In this way, the Dodecahedral conjecture is proved. 
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3 Computation 



The proof of the Dodecahedral conjecture is based on a series of computer calculations. 
This section briefly describes the computer algorithms, the code implementing those 
algorithms, and issues of the reliability of the computer code. 

There are three main computer programs that are used in the proof. The first is a 
graph generator that generates, up to isomorphism, all planar graphs satisfying a list of 
properties. The second is a linear programming package. The third is a piece of code 
based on interval arithmetic that automatically proves nonlinear inequalities over the 
real numbers. This section discusses each in turn. 

This section also discusses some additional computer programs. Although these 
computer programs, strictly speaking, are not part of the proof, they are relevant to 
understanding the structure of the proof and the reliability of the computer implemen- 
tation. We include a brief discussion of nonlinear optimization software, Tarski's deci- 
sion procedure for real-closed fields, and formal theorem proving packages. 

3.1 Electronic resources 

A permanent archive has been set up for all of the external resources related to this 
proof (26]. This archive is under version control by Google Code HI OH . The site consists 
of a download area where one may obtain the source code and supporting documents 
to this paper. Additionally there is a subversion [8| repository. This means that the 
snapshot of the code and documents in the exact form they took at the time of creation 
is permanently available. It also means that any changes (for instance, a bug fix) will 
leave a permanent public electronic trail. The major components of ll26ll are: 

1. The source code for the three different programs used in the proof (graph gener- 
ation, linear programs, and interval arithmetic inequality prover) 

2. A list of all tame Voronoi graphs, up to isomorphism. 

3. A list of inequalities that have been established by interval arithmetic. 

4. A list of the inequalities that have been used in linear programming. 
Additional information, such as source code documentation, is available as well. 

3.2 Tarski arithmetic 

The Dodecahedral conjecture, after a few preliminary reductions, can be expressed as a 
statement in the elementary language of the real numbers. The elementary language of 
real numbers is a first order language built from quantifiers (V, 3), logical connectives 
(A, V, =>, -t), functions symbols for ring operations (+, x, -), variables Xi, and constant 
symbols (0, 1). By a fundamental result of Tarski, the elementary theory of the real 
numbers is decidable. Thus, the truth of the Dodecahedral conjecture can be decided in 
theory by standard algorithms such as Collin's cylindrical algebraic decomposition Q, 
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or the Cohen-Hormander algorithm |2T|. However, in practice, these decision proce- 
dures take exponential time in the number of quantifiers, and thus are far too slow to 
be of practical value for this conjecture. 

To formulate the Dodecahedral conjecture as a statement in the elementary theory 
of the reals, consider Voronoi cells without truncation, centered at the origin. The vol- 
ume of the regular dodecahedron is an algebraic number u>dod, hence definable in Tarski 
arithmetic. Also, there are a priori bounds n on the number of faces of the Voronoi cell, 
and thus also on the size of the clusters of balls that give candidates for counterexam- 
ples. For example, by adding balls to the packing to decrease the volume, the packing 
becomes saturated. Assuming saturation, E. Harshbarger gives a quick calculation of 
n < 58 faces Ifl9l . (The exact value of this constant is not important as long as it is 
explicitly given.) The assertion of the Dodecahedral conjecture is then expressed as an 
enormous conjunction of cases, with conjuncts indexed by an explicit enumeration of 
all possible combinatorial structures of a Voronoi cell, including a fixed triangulation of 
each face of the cell. The fixed triangulation determines a partition of the Voronoi cell 
into tetrahedra. The volume of each tetrahedron is expressed by means of a Cayley- 
Menger determinant as a definable function of the the edge lengths of the tetrahedron. 
Thus, for each combinatorial structure X with m < 58 faces, a Tarski statement asserts 
that for all vectors vi, . . . , v m e R 3 , if the Voronoi cell defined by the vectors v, has 
combinatorial structure X, then the volume of the cell is at least tOdod- The outer block 
of 3 • 58 universal quantifiers - not to mention the nested existential quantifiers - is 
hopelessly beyond the practical reach of current algorithms. The statement that the 
regular dodecahedron is the unique minimizer can be similarly expressed. 

If the entire Dodecahedral conjecture can be expressed in Tarski arithmetic, then 
perhaps it is not so surprising that many of the intermediate steps in the proof can 
also be so expressed. These intermediate steps also tend to be beyond the reach of 
current decision procedures. But here the situation is not so hopeless. Many of these 
intermediate problems involve no more than a dozen quantifiers. One can imagine the 
day that these problems might fall within the reach of decision procedures for Tarski 
arithmetic. 

Describing various intermediate steps of the proof as exercises in Tarski arithmetic 
is a useful point of view. Doing so identifies a family of subproblems that can be 
expressed in a common language, and that can often be solved by similar techniques. 
The complexity of the problems can be measured objectively by counting the number 
quantifiers. 

Here are some geometrical objects that are definable within the Tarski arithmetic 
that appear in the proof of the Dodecahedral conjecture. By Heron's formula, the 
circumradius rj(x, y, z) of a triangle with sides x, y, z is elementary definable. Set 



for three points u, v, w e R 3 . The volume of a tetrahedron with vertices Vi, . . . , V4 is 
elementary definable. In fact, by Cayley-Menger determinants, there is a polynomial 
A of the six squared edge lengths such that 



T] V (u,v,w) = rjQu - v|, |v - w\, \u - w|), 
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In particular, A(xy) > 0, whenever the variables Xy have the form |v, - \>j\ 2 for some 
points Vj e R 3 02J Lemma 8.1.4]. 

The orientation of V\ in S — [v\, x>2, V3, V4} is said to be positive (zero, negative) 
when Vi and the circumcenter of S lie on the same side of (resp. on, on opposites 
side of) the plane passing through {V2, V3, V4} [ 18, Lemma 5. 15]. This is an elementary 
condition. 

For S — { vi , . . . , v,-} c R 3 , let aff + (0, S ) be the cone with apex generated by S : 

{fiVi +■•• + t r v r : ti > 0). 

Write aff^OjS) for the corresponding set with strict inequality f, > 0. Let conv(S) be 
the convex hull 

r 

{hvi +---+t r v r : ti > 0, Yj'i = !)■ 

!=1 

Write conv°(S ) when the inequality is strict f, > 0. 
Define the right circular cone rcone by the formula 

rcone(v, w, h) — \x : (x — v) • (yv — v) > |jc — v| \w - v\ h}. 

For a set S of fixed finite cardinality, the membership conditions x e aff + (0, S), x e 
cone(S), x e rcone(S) can be expressed in the Tarski language. 

In preparing this abridged version, the proofs of numerous statements in the Tarski 
language (involving a small number of quantifiers) have been omitted. These tend to 
be the arguments that can most easily be skipped without disrupting the overall flow of 
the proof. The full proofs appear at lfT6l . 

3.3 Formal proof 

A formal proof is a proof in which every logical inference has been checked, back to 
the foundational axioms of mathematics. Except in trivial cases, a computer is used 
to generate a formal proof, because of the large number of inferences involved. Both 
conventional proofs and computer assisted proofs can be formalized. In a computer 
assisted proof, this amounts to a formal verification of the correctness of the computer 
code. Formal verification of computer code is a difficult task. For that reason, formal 
verification tends to be reserved for situations where correct performance is critically 
important, such as the verification of aircraft control systems, cryptography algorithms, 
and security protocols for the internet. There is no other means of checking computer 
software that can assure reliability at levels that remotely compare with the assurance 
afforded by formal verification. 

There is a long term project, called Flyspeck Light, to give a formal verification 
of the Dodecahedral conjecture. Although this project is far from complete, parts of 
this project have already been carried out. This means that some of the computer code 
for this project now carries a proof of correctness, according to formal mathematical 
standards. One such program is discussed in the next subsection. 
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3.4 Graph generation 



The proof of the Dodecahedral conjecture is based on three separate computer pro- 
grams. The first of these is a planar graph generator. It generates all planar graphs, 
up to isomorphism, that satisfy a given list of restrictions. The restrictions include a 
bound on the number of vertices in the graph, so that it is obvious that only finitely 
many graphs are possible. 

The correctness of the graph generator program has been the subject of extensive 
mathematical investigation. Early versions of the program were written by Hales in 
1994 (in Mathematica), 1997 (in Java), and 2000 (in Java). The same computer pro- 
gram is used in the Kepler conjecture and Dodecahedral conjecture. They differ only 
in their input parameters. This computer program became the subject of G. Bauer's 
dissertation in computer science at the Technical University of Munich |2|. This 172- 
page dissertation translates the Java code into the formal theorem proving system Is- 
abelle/HOL Oil and gives a detailed mathematical treatment of the graph theory un- 
derpinning the computer code. The dissertation analyzes every line of code. Building 
on the work of this thesis, B. Bauer and T. Nipkow have completed the formal correct- 
ness proof of the HOL implementation of the graph generator ||29l . (Their published 
article mentions only the Kepler conjecture, but the formal verification has been ex- 
tended to apply to the input parameters of the Dodecahedral conjecture as well.) As 
a result of this work, the graph generation is currently the most scrupulously checked 
part of the proof of the Dodecahedral conjecture. 

There are several published sources that provide details of this algorithm, and there 
is no need to repeat details here. The basic idea is to start with a small set of planar 
graphs (called 'seed' graphs) with the property that every planar graph to be classified 
is known to have one of the seed graphs as a subgraph. The seed graphs are then 
extended by adding one face at a time. Faces are added in all possible ways so that 
it is clear at every step of the algorithm that every biconnected planar graph will be 
generated. At the same time, pruning operations discard partially completed graphs 
when it can be shown that the partial graph is not a subgraph of any of the graphs to 
be classified. The pruning operations prevent a combinatorial explosion of cases. See 

ma §5], ma §19], a usa. 

3.5 Linear programs 

The second computer program that is used in the proof of the Dodecahedral conjecture 
is linear programming. There are several hundred linear programs that appear in the 
proof. 

Formal verification has not yet been extended to this portion of the computer code. 
However, the recent dissertation of S. Obua takes the first steps in this direction 11501 . 
That work gives a formal correctness proof of the basic linear programs that appear in 
the proof of the Kepler conjecture. In particular, he has developed all the infrastructure 
needed to carry out formal correctness proofs of linear programming problems. The 
formalization completed by Obua is a larger project than the formalization of the linear 
programming segment of the Dodecahedral conjecture. Thus, one can expect that the 
formalization of this piece of computer code will soon follow suit. 
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In the 1998 proof, computer code written in C generated the linear programs, which 
were then fed to the commercial linear programming package CPLEX. In preparation 
for a formal proof, the computer code has been rewritten in the programming language 
Standard ML (SML), with an external interface to the solver GLPK. 

It is not necessary to trust the algorithms of the linear programming packages (such 
as CPLEX and GLPK) that solve the linear programs. These packages produce dual 
certificates that can be used to give independent verification of the solutions of the 
linear programming problems. 

In the proof of the Dodecahedral conjecture, the following situation arises. The 
objective is to prove that the maximum of a linear function x t-> cx is less than a given 
constant M, when subject to a system of linear constraints Ax < b. Here x, c, b are 
vectors with real entries and A is a matrix with real entries. The products are given 
by matrix multiplication of compatibly sized matrices and vectors. A vector inequality 
u < v means that w, < v, for each coordinate ;'. Explicit lower and upper bounds on the 
variables are given: ( < x < u. Expressed equivalently, the objective is to show that the 
linear system of inequalities 

Ax < b, I < x < u, cx> M 

has no solutions in x. The external linear programming package produces a dual cer- 
tificate in the form of a vector y, which that package claims to have the properties 

yA=c, y>0, yb < M. (1) 

If y indeed has these properties, then for any x satisfying Ax < b, it follows that 

cx = yAx <yb < M (2) 

as desired. 

Because of inexact arithmetic used by the external packages, these identities will 
only be approximately correct. The imprecision in the dual certificate can be readily 
eliminated as follows. If u is any vector, let u + be the vector obtained by replacing the 
negative entries of u with 0, and let uT be the vector obtained by replacing the positive 
entries of u with 0. By replacing the vector y with y + , the vector y satisfies < y. In 
the following lemma, e\ and 62 are small error terms that result from machine approxi- 
mation. By including them in the bounds on cx, a rigorous bound can be recovered. 

Lemma 3.5.1. Suppose that the real-valued vectors and matrices A,A\,A2, c,c\,C2, 
x, b, {, u, y satisfy the following relations 

Ax < b, A\ < A < A 2 , ci < c < c 2 , I < x < u, < y. 

Define residuals 

e\ = ci -yA 2 , e 2 = c 2 -yAi. 

V 

yb + ej M+ + e \ u + e 2 ' + + e i ' < M ' 

then cx < M. 



12 



Proof. S. Obua has given a formal proof of this lemma in the Isabelle/HOL system l30l 
3.7.2]. In fact, the proof follows from a simple embellishment of Inequality [2] 

cx — yb < c 2 x + + c\x~ — yAx — (e2X + + e\x~) + y(A\ — A)x + + y(A2 — A)x~ < ■ ■ ■ . 

□ 

The numerical data A\, A2, c\,C2, (, u,y, b are all explicitly given, so that the method 
yields explicit bounds. It is not necessary to trust the package that produces the certifi- 
cate y, because all of the assumptions of the lemma can be checked directly with simple 
matrix multiplications. The reliability of these matrix multiplications is guaranteed by 
using interval arithmetic. 

3.6 Interval arithmetic 

The third major computer program that is used in the proof of the Dodecahedral con- 
jecture is an nonlinear-inequality prover over the real numbers based on interval arith- 
metic. This subsection describes the methods involved and the computer implementa- 
tion. 

A finite number of nonlinear functions f\,..., f r are given. It is assumed that all 
functions have the same domain 

R = [auh] x [a 2 ,b 2 ] x [a m ,b m ], 

given as a product of intervals in R™, for some m. The computer program verifies that 

(/i (jc) > 0) v (f 2 (x) > 0) v • ■ • v (f r (x) > 0), (3) 

for every point x e R. The approach is similar to the approach described in R. B. 
Kearfott [24], based on interval arithmetic. Our methods are similar to algorithms in 
widespread use for rigorous global optimization. Closely related algorithms are also 
described in 11371 . 

The method is based on an iteratively refined partition of the domain R into a finite 
number of smaller and smaller rectangles that cover R. 

Start with X = {R}, then repeat the following procedure. Pick T e X; replace X 
with X \ {T}; and calculate a lower bound fj(x) > cij(T), for all x 6 T. If aj{T) > 
for some /', then the desired bound Q holds on T. If the desired bound holds on T, 
then continue to the next rectangle in X. Otherwise, choose finitely many rectangles, 
7i , . . . , T/t that cover T; replace X with X U {Ti, . . . , 7^}, then repeat. The procedure 
terminates when X — 0. (If the procedure is applied to a false disjunction <j3j>, there is 
no termination.) 

When subdividing rectangles to obtain smaller covers, one does it in such a way the 
the width of the rectangle tends to zero as iteration continues. In this way, the lower 
approximations aj(T) to fj can be arranged to converge to the true minimum of fj on 
the rectangle. 

The lower bounds aj(T) to a function fj on a rectangle T are obtained by methods 
of interval arithmetic. The function / = fj is generally C 2 on its domain|^The function 

2 The functions encountered in practice are usually C 2 , but not always so. When functions are not C 2 , 
Taylor approximations are avoided. 
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/ can be expanded in a Taylor polynomial approximation with explicit error bounds. 
Derivatives are calculated by automatic differentiation. The error bounds are based on 
the Lagrange form of the error term in the Taylor approximation. Interval arithmetic 
produces rigorous bounds on the error terms. 

There have been several separate implementations of the interval arithmetic pack- 
age. The source code for all of these packages is publicly available. (The code that 
is used for the proof of the Kepler conjecture is the same as the code that is used for 
the proof of the Dodecahedral conjecture. Only the statements of the inequalities to 
be proved differ.) The first version, written in C++, was developed by T. Hales over 
the period 1994-1998. A second version, written in C, was developed by S. Fergu- 
son 1995-1997. A third version, written in SML, was developed by S. McLaughlin 
2006-2008. A interval arithmetic package has also been developed by R. Zumkeller 
for the theorem proving system Coq [3], although it has not been used to give a formal 
verification of any of the inequalities that arise in the proof of the Dodecahedral con- 
jecture ||38l . These implementations are all independent of one another. (Algorithms 
were shared among us, but the code was independently implemented.) By comparing 
the proofs of different inequalities in different systems, we have developed a high de- 
gree of confidence that the implementations of the algorithms are essentially correct. 
Of course, it would be desirable to have a formal correctness proof, but this part of 
the Flyspeck Light project has not been completed. (The SML implementation and 
Zumkeller's research are partial steps in this direction.) 

The list of nonlinear inequalities that are used in the proof of the Dodecahedral con- 
jecture appears at |26 1. The domains of the functions are subsets of R m , for m < 7. The 
complexity of the verification increases rapidly with m. This proof implements some of 
the tricks introduced in [ 18] to reduce the dimension of the domain wherever possible. 
Dimension reduction is based on established monotonicity properties of the functions. 
(For example, the volume of a Voronoi cell does not increase when it is intersected with 
a half-space.) Whenever the functions fj are twice continuously differentiable, a first 
order Taylor polynomial with explicit error bounds on the second derivatives is used. 

The functions fj represent elementary geometric quantities such as linear combina- 
tions of angles, dihedral angles, solid angles, and volumes. Explicit formulas for these 
functions are known involving rational functions, the square root, and arctan functions. 
The typical form of a function fj is a linear combination of terms of the form 

arctan(a/ Vfr), 

where a, b are explicit polynomials on R m . In the formula d3j, the number of disjuncts 
is usually one r = 1, but in some cases there are two disjuncts. 

The computer calculations use interval arithmetic to control for floating-point round- 
ing errors. Every real number x is represented on the computer as an interval [a, b] con- 
taining x, where a and b are exactly representable floating point numbers Q][32). The 
calculations conform to IEEE-754 standards [23]. Approximations to inverse trigono- 
metric functions are based on published approximations GUI . 
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3.7 Nonlinear optimization 

Previous subsections describe the three main pieces of computer code used in the proof 
of the Dodecahedral conjecture: graph generation, linear programming, and interval- 
arithmetic inequality proving. This subsection describes one additional software pack- 
age indirectly involved in the proof: nonlinear-optimization. 

The disjunction of inequalities in formula ([3]) can be represented as a constrained 
minimization problem: show that the global minimum of f\ on the domain 

{x e [a u bi] X ■ ■ • X [a m , b m ] : f 2 (x) < 0, . . . , f r (x) < 0} 

is positive. 

Nonlinear optimization libraries have been used to test all the inequalities in the 
collection ||25l . |6). The code generates a large random set X of points in the domain 
and runs the algorithm for each initial point x e X to find a local minimum to the ob- 
jective function f . If X is sufficiently large and sufficiently random, it can be expected 
that one of the local minima produced to be a numerical approximation of the global 
minimum. 

In practice, this approach works remarkably well on this collection of problems, 
largely because the functions fj tend to be rather bland from the point of view of non- 
linear optimization. (Typically, the second derivatives of fj are small; the level surfaces 
of fj are approximately planar; there are no local minima in the interior of the domain; 
the global minimum occurs at a corner of the domain; and every run of the algorithm 
produces the same local minimum.) Thus, the method can usually determine the true 
global minimum with high probability. 

If this nonlinear optimization is not part of the proof tree, what purpose does it 
serve? First of all, although we have tried to be careful to avoid any errors in the 
computer code, an independent check of the results is certainly welcome. It makes the 
proof more robust against possible errors. In fact, this independent check has helped 
us to spot and correct data entry errors. Secondly, the package was used to discover 
inequalities that were likely to be true, and to discard quickly inequalities that were 
false. The plausibly true inequalities then became candidates for rigorous nonlinear 
optimization with interval arithmetic. 

3.8 Summary 

This section has described various computer programs and algorithms that have been 
used in the proof of the Dodecahedral conjecture. The rest of the article assumes that 
the types of computations described in this section can be reliably performed by com- 
puter. 

Recall briefly how the three main computer programs enter into the proof. A planar 
graph is associated with each potential counterexample of the Dodecahedral conjecture. 
The properties of this graph is studied, and it is shown to be a tame Voronoi graph. Us- 
ing the graph generator program, all such graphs are classified up to isomorphism. This 
reduces the proof to a finite enumeration of cases. Linear programs are then used to 
show that each case in this enumeration is infeasible. The nonlinear inequalities appear 
in several different parts of the proof. They are used, for example to establish that the 
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graph associated with a counterexample is a tame Voronoi graph. Nonlinear inequal- 
ities are also used to justify the list of inequalities used in the linear programs. (The 
linear programming inequalities come as linear relaxations of nonlinear inequalities.) 

4 Fejes Toth's Reduction 

L. Fejes Tdth proved the Dodecahedral conjecture under the extra hypothesis that 
A(vo, ltdod) \ ( v o) has at most 12 elements. The proof occupies about eight pages of the 
book Regular Figures 1351 . See the main theorem of Section 33 and the main theorem 
of Section 41 (including the note on page 265). Fejes Toth's bound is a general bound 
about the volume of truncated polyhedra. The polyhedra do not need to be Voronoi 
cells in a sphere packing. Here is a sketch of his proof. 

Theorem 4.1. Let P be any polyhedron with at most 12 faces that contains a unit 
sphere S 2 . Let B be the ball of radius tdod concentric with S 2 . The volume of the 
intersection of P with B is at least u>dod- Equality holds exactly P is congruent to the 
regular dodecahedron Q.(Ad d, 0). 

Proof. (Sketch) By translation, the proof reduces to the case that the origin is the center 
of S 2 . If k < 12, there is a polyhedron with smaller volume and k + 1 faces, obtained 
by clipping a corner of the polyhedron with a new face. Thus, assume k = 12. View 
the case in which some vertices of the polyhedron have degree greater than 3 as degen- 
erate cases of polyhedra where all degrees are three, where some of the vertices have 
coalesced. With these conventions, there are 12 faces, 30 edges, and 20 vertices. 
Let 

. . _ 1 5 sec 2 (*), x < 6 

8X ~{±sec 2 (0 o ), x>6 , 

where 6q is defined by sec(#o) = td d- For each face F, of the polyhedron, let w,- be 
the point on S 2 closest to the plane through F,. Let Si be the radial projection of Fj to 
a spherical polygon on S 2 . The volume a> of the fj rf-truncated polyhedron satisfies 

12 „ 

w>J] J g(8(w h x))dx (4) 

!=1 

where dx is the usual measure on S 2 , and 6{wi,x) is the geodesic length of the arc 
on S 2 joining w,- to x. The integral on the right is exactly the volume of the truncated 
polyhedron obtained by projecting each polygon S, back out to a plane through w t 
parallel to Fj. In particular, equality holds if the plane through each face is tangent to 
S 2 . 

By the estimate of Section 33 of the book, the integral on the right is at least 

120 j g(0(w,x))dx, 

where T is a spherical triangle with angles n/2, 7r/5, 7r/3. Here w is the vertex of 
the triangle T that has angle n/5. (This estimate holds more generally for any non- 
decreasing function g.) 
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This integral is precisely the volume of a regular dodecahedron of inradius 1. In- 
deed, when the polyhedron is a regular dodecahedron, the maximum of 8(w, x) over 
the pentagon Si is exactly 6>o. (This fact is equivalent to the definition of tdod as the cir- 
cumradius of the regular dodecahedron, so that no truncation occurs.) So g(\w - x\) = 
(sec 2 0(w, x))/3. The inequality is an equality. Also, each regular pentagon 5, can 
be triangulated into 10 triangles with angles 7t/5, n/2, x/3. The 12 faces then give 120 
triangles T. The result follows. 

L. Fejes Toth also considers the case of equality and finds by similar arguments 
that the only minimizing polyhedron is the regular dodecahedron. □ 

5 Geometry of Voronoi Cells 

This section describes the basic geometry of the Voronoi cell and its truncation. 

5.1 Basic truncation 

Let A = A(0, Itdod) be a finite packing containing 0, and let f2(A, 0) be the Voronoi 
cell. Let B(x,r) be the closed ball of radius r centered at x e R 3 . Let £lo(A, 0) = 
Q(A, 0) n B(Q, t,j ot i). Also, Q fra „ c (A, 0) has been defined in Section [2~T| 

Recall A* = A \ {0}. There is a graph G(A) with vertex set A*, whose edges 
are formed by pairs [v,w] such that < |v — w| < ltdod- Figure [3] shows the regular 
icosahedron, which is the graph G{A c f„d)- 

5.1.1 One cap 

The geometry of Oo(A, 0) will be discussed first, then adapted to fi, rM „ c (A, 0). The 
truncated cell Oo(A, 0) is obtained from the ball B(0, td d) by removing a spherical cap 

C(v, tdod) = {x e B(0, t dod ) : \x - v| < \x\] 
for each v e A*. Each spherical cap is bounded by a sphere of radius tdod an d a planar 
disk formed by the intersection of the bisector of {0, v) with the ball B(Q, tdod)- 
Volume and solid angle calculations use the following functions: 

4>(h, t, A) = A v th(t + h)/6 + A s 

A(h, t,A) = (1 - h/t)((/>(h, t, A) - 4>(t, t, A)), ' 

where A = (A v , A s ). The subscripts v and s on the components of A stand for 'volume' 
and 'solid angle.' This terminology is justified by the following calculation. When 
A* = {v} has a single point, the truncated Voronoi cell is a ball of radius tdod with a 
single cap C(v, tdod) removed. Its volume depends only on |v] = 2h (and tdod)- The cone 
at the origin generated by the points of C(v, tdod) is rcone(v, hjtdod)- An elementary 
calculation gives the volume and solid angle formulas. 

vol(C(v, tdod)) = ^(h,tdod,(l,0))so\(icone(v,h/tdod)) (f . 
<p(h,t d od,(0,D) =1. } 

Thus, <p{h, t, A) is the ratio converting arbitrary linear combinations of volume and solid 
angle into solid angle. 
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Figure 4: Quoins 



5.1.2 Two caps 

When A* = {v, w) contains two points, the two spherical caps meet if and only if 
j] V (Q,v,w) < tj oc i. (Recall that j] V denotes the circumradius.) The volume formula 
for Oo(A, 0) as a function of v,w is continuous across the hypersurface 77,(0, v,w) = 
t dod , but not analytic. When r] v (Q, v, w) > tdod the caps are disjoint and the volume is 
independent of |v - w\, depending only on |v] and |w|. That is, the volume does not 
depend on the location of the caps, provided they are disjoint. 

When rj V (Q, v, w) < tj j, the volume depends on |v|, |w|, \v—w\. Note that rj v (Q, v, w) < 
tdod implies |v - w\ < 2t dod , so that the graph G(A*) contains the edge {v, w}. By 
inclusion-exclusion, the volume of Q.q(A, 0) is 

vol(fi(0, t dod )) - vol(C(v, t dod )) - vol(C(w, t dod )) + vol(C(v, t dod ) n C(w, t dod )). (7) 

The set C(v, t dod ) n C(w, t dod ) can be partitioned into four regions (called quoins). 

Let Pi = aff{0, v, w] be the plane through 0, v, w. Let P2 be the plane orthogonal to 
Pi that passes through and the circumcenter of the triangle {0, v, w). Then R 3 \ (P t U 
P2) contains four connected components, partitioning C(v, tdod) n C(w,t dod ) into four 
quoins. Let a = |v|/2, b = rjv(0, v, w), c = tdod- Let q(a, b, c) be the volume of the quoin 
given by the intersection 

C(v, tdod) n C(w, frf orf ) r\H x ^H 2 , 

where //] is a half-space bounded by P\ and H2 is the half-space bounded by P2 con- 
taining w. The volume of a quoin is computed in 1 18 §7.3]. If a < b < c, then 

6 quo(a, b, c) — (a + 2c)(c - a) 2 arctan(e) + a(b 2 - a 2 )e 

-4c 3 arctan(e(/3 - a)/(b + c)), where (8) 
e\b 2 - a 2 ) = (c 2 - b 2 ), and e > 0. 

Otherwise, quo(a, b, c) = 0. By Q, this gives the volume of ^o(A, 0) when A* = {v, w}. 



5.1.3 Three caps 



Lemma 5.3.3 shows it is not possible for four or more spherical caps to meet. Here we 
consider the case of three spherical caps. 
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When A* = {vi,Vz, V3} contains three points, the three spherical caps meet if and 
only if the circumradius of the simplex S = {0, Vi, Vz, ^3} is less than tdod- When this 
happens, each edge satisfies |v; - Vj\ < 2td d, so that the graph G(A) is a triangle. Also, 
the circumradius of each face is at most tdod < V2. Thus, 

{0,vi,V2,v 3 }e5(A,0). 

The convex hull conv(S) is singled out for special truncation in Cl, nmc . 

If the three spherical caps do not meet or meet in at most pairs, then the volume, 
by inclusion-exclusion is given by a formula similar to (|7]). In particular, the volume is 
given in terms of quoins, and so forth. 

The case when the graph G(A) is connected, but not a triangle has particular inter- 
est. Suppose that there is no edge between Vi and V3. In this case, the volume of £2 ( ) 
depends on |v,-|, |vi - V2I and |v2 - V3 1, but not on |vi - V3I. Thus, the point V3 can be 
moved, subject to the constraints fixing |v3 1 and |v2 - V3I without changing the volume. 
In particular, V3 can be moved until | V3 — v 1 1 = 2td„d- This results in the following simple 
lemma. 

Lemma 5.1.1. Suppose G(A) contains three vertices, is connected, but not a triangle. 
Then there is another packing A' whose graph G(A') is a triangle and such that 

vol(Q (A,0)) = vol(Q (A',0)). 

The same lemma holds for Q. mmc because in the situation at hand Qq = Qtrunc- 

5.2 Planarity 

Definition 5.2.1 (spherical). For each edge e — {v,w} of a graph G' on a vertex set 
A*, let A e be the arc on the unit sphere at the origin formed by the intersection of the 
sphere with aff^O, {v, w}). The graph G' is said to be spherical if the sets A e do not 
meet one another, as e runs over the edges ofG'. 

A spherical graph is clearly a planar. 

Lemma 5.2.2. G(A) is a spherical graph. 

Proof. It suffices to show that the sets aff^(0, {v, w}) do not meet one another, as e — 
{v, w] runs over edges. This is a statement in Tarski arithmetic. A detailed proof ap- 
pears in |[T6l Lemma 3.2]. It is the reparametrization of a published theorem lfT3l 
Lemma 3.10]. □ 

5.3 Triangles 

This subsection describes the geometry associated with triangles in the graph G(A). 
For each triangle, there is a set {0, vi, V2, V3} such that 

hi < 2w, \Vj - vj\ < 2t dod , i, j < 3. (9) 
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Lemma 5.3.1. Let S = {0, Vi , V%, V3 } be a set of four points that satisfies ([9]). Then there 
does not exist w € conv(5 ) that satisfies \w — u\ > 2 for allu 6 S. 

Proof. This is a statement in Tarski arithmetic that can be expressed with 12 quantifiers 
(3 coordinates for each of Vi, Vi, V3, w). The proof is the reparametrization of lfT8l 
Lemma 4.15]. It is written out in full in 1 16, Lemma 3.3]. □ 

Lemma 5.3.2. Let S = {0, v; , V2, vj} be a set of four points that satisfies ([£]). Then there 
does not exist w € aff + (0, {vi , Vi, V3}) ( the cone with apex generated by positive linear 
combinations ofv t ) that satisfies \w — u\ > 2 for u e S and \w\ < 2tdod- 

Proof. This is a statement in Tarski arithmetic. The case when w e conv(S) is covered 
by the previous lemma. The remaining case is a reparametrization of lfT8l Lemma 4.19]. 
See also, |[l6l Cor 3.7]. □ 

Lemma 5.3.3. Let 0, v\ , v^, V3, v\ e A, with v,- distinct. The intersection nf =1 C(v,-, tj oc i) 
of spherical caps is empty. 

Proof. A nonempty intersection implies that each edge {v,, vf] belongs to G(A), form- 



ing a complete graph on four vertices. The graph is planar by Lemma 5.2.2 Its planar 
representation is a triangle with one vertex inside, connected to all three vertices of the 
triangle. Geometrically, this corresponds to a point V4 e aff + (0, {vi, V2, V3}), which is 
impossible by the previous lemma. □ 

Lemma 5.3.4. Let [v\,V2, V3} and (v'pV^vU be two distinct triangles in G(A). Then 
aff^O, {vi, V2, V3}) is disjoint from affl(0, {v' p v' 2 , v^}). In particular, conv({0, Vi, V2, V3}) 
is disjoint from conv({0, v' p vi, v' 3 }). 

Proof. Since the graph is planar, if the two cones intersect, then one triangle must be 
contained in the other triangle. That is, one cone is contained in the other. This leads 
to a vertex w of one triangle in the other cone. This is prohibited by Lemma 5.3.2 □ 



As mentioned earlier, a three-fold intersection of spherical caps produces a triangle 
in the graph G(A) and a set 

{0,vi,v 2 ,v 3 } eS(A,0). 
The next two lemmas investigate the geometry of such {0, Vi , Vz, V3}. 

Lemma 5.3.5. If S — {0,Vi,V2,V3} is a set of four points such that each face has 
circumradius at most V2 and such that \u—v\ < 2td oc ifor u,v € S , then the circumcenter 
of S lies in conv(S). Also, if x e conv(S) and if some w has distance at least 2 from 
each point of S , then x is at least as close to some point of S as to w. 

Proof. This is a statement in Tarski arithmetic. The statement about the circumcenter 
is lfT8l Lemma 5.18]. (No reparametrization is needed.) If x comes closer to to w than 
to each point of S , then the Voronoi cell Q.(S U {w}, w) meets conv(S). Again by [18, 
Lemma 5.18], this implies that the circumradius of some face of S is greater than V2, 
which is contrary to hypothesis. An alternative proof of both parts of the lemma is 
contained in ifTBI Lemma 3.5.3.61. □ 
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The following lemma justifies limiting packings to those satisfying A = A(0, Itdoi)- 
Extending the packing A beyond radius 2tdod cannot decrease the volume of the trun- 
cated Voronoi cell. 

Lemma 5.3.6. Let A = A(0, 2t dod ). Let A c A', where A'(0, 2t dod ) = A(0, 2t dod ). Then 

O (ra „ c (A',0) = Q 

I rune 

(A,0). 

That is, the truncated Voronoi cell cannot be decreased in volume by adding additional 
points to the packing outside the ball B(0, 2t dod ). 

Proof. Let w e A' \ A (so that |w| > 2t dod ). Every point in B(0, t dod ) n Q, trunc (A, 0) is 
clearly closer to than to w so belongs also to Q,„ me (A', 0). 

Assume that x t B(Q, t dod ) and x e Q fr „„ c (A, 0). Then by the definition of Q fnmc , 
there is some S = {0, Vi,v 2 , V3} e S{A,Q) such that x e Q(A,0) n conv(S). By 
Lemma [5. 3. 5| x is closer to than to w. Thus, x e Q.(A',Q). The result follows. □ 



Assume S = {0, Vi, V 2 , V3} e S(A, 0). From Lemma 5.3.5 it follows that 



aff° (0, {vi , V 2 , v 3 }) n Q(A, 0) = aff^O, {vi, v 2 , v 3 }) n Q /ra „ c (A, 0) = conv(S ) n Q(S, 0). 

(10) 

That is, the calculation of volume can be made locally in S without reference to the po- 
sition of the packing A. A formula for this volume calculation appears in [13, §8.6.3]. 
This gives the formula in the case when the circumradius is at most t dod . When the 
circumradius is at least t dod , the spherical caps intersect in pairs, and the inclusion- 
exclusion formula (using quoins) can be used |7]). In summary, the inclusion-exclusion 
formula can be used for all calculations of truncated Voronoi cells Q. tnmc (A, 0), except 
for {0, Vi, v 2 , V3} e «S(A,0) with a circumradius less than t dod . In this case, the explicit 
formula just mentioned applies. 



5.4 Connecting the graph 

The volume formula for Q. tnmc (A, 0) depends only on |v|, for v e A* and on |v - w| for 
{v, w) and edge of G(A). If the graph G(A) is not connected, the collections of spherical 
caps for two different connected components of the graph do not intersect one another. 
Thus they form two or more non-interacting "islands" of spherical caps that can be 
moved independently around the globe B(0, t dod ) without changing the volume of the 
truncated Voronoi cell. In particular, one island of spherical caps can be moved rigidly 
until some vertex v of on connected component of the graph has distance exactly 2t dod 
from some vertex in another component. This connects two components of the graph 
without changing volume. Thus, every truncated Voronoi cell has the same volume of 
another with a connected graph. Assume without loss of generality that the graph is 
connected. 

A biconnected graph is defined as a connected graph that has no articulation ver- 
tices; that is, by removing any vertex, the graph remains connected. Consider a con- 
nected graph G(A) that is not biconnected. Then there exists some w e A* that dis- 
connects the graph. Write the vertex set A* as a disjoint union of three sets {w}, Aj, 
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A2 such that there is no edge between A; and A2. The volume of 0, nmr (A, 0) is inde- 
pendent of the distances between points of A] and A2. This means that A2 can move 
rigidly, while constrained to preserve |v] for v e A2, |v - u\ for u, v e A2 U {w}. Continue 
the rigid motion of A2 until some distance \u - v\ decreases to 2fj ^, for u e A] and 
v e A2. (See Lemma 5.1.1 ) Repeat this construction until G(A) becomes biconnected. 
The process does not alter the volume of Q fnmc (A, 0). Now assume, without loss of 
generality, that G(A) is biconnected. 



5.5 Standard components 

This subsection continues to assume the standing list of assumptions on A. Specifically, 
e A = A(0, tdod), and G(A) is a biconnected planar graph. Under these assumptions, 
each face of G(A) is a simple polygon. In particular, there are no vertices of degree one 
in the graph. 

For the moment, the situation can be generalized somewhat to allow G' to be any 
graph on the vertex set A* that is spherical and biconnected. For each edge {u, v} of G', 
form the cone aff°(0, {u, v}). Let X = X(G') be the union of these cones and let T(G') 
be the complement of X in R 3 . The open set F(G') breaks into a finite set of connected 
components. Write [T(G')] for this set of connected components. The set [T(G')] is 
in natural bijection with the set of faces of the graph G'. Write Uf for the connected 
component corresponding to face F. The component U is said to be indexed by F. 

Each face F is identified with a sequence of vertices (vi,...,v r ) with v, e G', 
giving the cyclic order of the vertices around the face. The sequence is well-defined up 
to cyclic permutation, so that (vi, . . . , v r ) defines the same face as (v,., Vi, . . .). Pick the 
order of the cycle counterclockwise around each face. 

Each connected component U has a solid angle sol(f/), which is defined to be the 
area of U n S 2 . The sum of the solid angles is the area of S 2 : 

so\(U) = 4-n. 

Ue[Y(C')] 

If U = Uf is a connected component and v a vertex of F, then there is an azimuth 
angle assigned to (U, v) with the property that if F\, . . . ,Fj are all the faces of G(A) 
that contain v, the sum of the azimuth angles around v is 2n: 

k 

^ azim(t//r., v) = 2n. 

(=i 

By definition, the azimuth angle equals the interior angle of the spherical polygon 
U C\ S 2 at v/\v\. By Girard's formula for the area of a triangle or polygon, for F = 
(vi,...,v r ): 

k 

so\(U F ) + (r-2)7i = Yj azim(f/ F , v ; ) (11) 

Set u>(A) = vol(Q fraflc (A, 0)). If U is a connected component of T(G'), set 
w(A, U) = vo\(U n Q fr „„ c (A, 0)). 
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Then 

w(A)= J] "(A, I/). (12) 

£/e[Y(G')J 

Set M ( ; 0( ; = 0.42755 and set //(A) = <w(A) - AnMdod- The desired inequality can be 
expressed as /i(A) > p.(Ad„d)- We have 

MA^) * 0.177540. 

The number pXAdod) is called the squander target in [16]. When £/ is a connected 
component of T(G'), set 

//(A,t/) = w(A,C7)-M dod sol([7). (13) 

Then 

/i(A) = J] MA, CO. 

£/£[7(A)] 

Now specialize again to the situation where G' = G(A). In this case, write X(A) = 
X(G(A)), Y(A) - Y(G(A)), and so forth. A connected component of F(A) is called a 
standard component. (This term is a reparametrization of a term by the same name in 
the proof of the Kepler conjecture.) 

If t/is indexed by a triangle F - {v\, V2, V3} in the graph, then U = aff^O, {vi, v%, V3}). 
In this case, cl>{A, U) is precisely the volume of the region already considered in Equa- 
tianQjj] 

If U is not indexed by a triangle, then 

w(A, U) = vol(t/ n Q (A, 0)). 

The formula for w(A, t/) in this case follows by inclusion-exclusion as in Section [5TT] 
Suppose that F is a face of G(A) whose vertices are given by (vi, . . . , v r ) (listed con- 
secutively around the face). Set h/ = |v,|/2, t = t^, b* = rj v (Q, v,-, v 1±1 ), A = (1,0). By 
Ifl8l Eqn. 7.12], inclusion-exclusion gives 

cj(A,Uf) =so\(U F )(f>(t,t,A)+ 

X - =1 (azim(£/F, v,)A(/i ; , t, A) + quo(/!,-, ^7+, t) + quo(/i,-, for, f)). 

The derivation of this formula relies on two geometric facts. First, each quoin lies 
entirely in a single standard component. Second, for each 1 < i < r, let P ± be the open 
half-plane bounded by the plane through {0, v,, v, ± i), so that P+ n P n V = Uf C\V, 
for some neighborhood V of v. Then, C(v,) n P+ n P_ c f/f . These facts are justified 
in lfT~8l Lemma 12.5]. The reparametrized version appears in |[T6l . 

Theorem |6. 1 1 shows that //(A, f/) is positive for every standard component U. The 
constant Mdod is chosen so that the minimum of /((A, U) - as both A and U vary - is 
very close to zero (about 10~ 7 ). The function fi tends to have better numerical behavior 
than at. For that reason, even though the two functions carry essentially the same 
information, estimates are expressed in terms of fi rather than a>, whenever possible. 
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Remark 5.5.1. The function p is closely related to a function t(-, t) that is used in the 
proof of the Kepler conjecture. The function p is, up to a small error term, a positive 
multiple of T(-,tdod)- The small error term comes from the fact that in this article, 
the constant Md od is used, and in [18] the constant Mq — l/(3(5, e/ ) is used, where 
6 tet — a/8 arctan( V2/5). (The constant 6 tet ~ 0.7796 is Rogers's famous bound on the 
density of sphere packings.) The difference is small: 

M -M dod ~ 1.86 x 1(T 7 . 

Because of the close similarity between p and t(-,?o), for every estimate involving 
r(-,?o) there is apt to be an analogous estimate involving p. The translation involves 
replacing Mq with Md c ,d, ?q with td„d and reseating the resulting function by an explicit 
positive scalar to get p. 

6 The Main Estimate 

This section proves the main estimate, which gives a lower bound on the function 
p(A, U) for any standard component U. The standing assumptions on A remain in ef- 
fect: € A = A(0, 2tdod) and G(A) is biconnected. This section makes no assumptions 
on the cardinality of A except where explicitly stated. 

Theorem 6.1. Let A. be a finite packing satisfying the standing assumptions. Let U p 
be a standard component indexed by a face F ofG(A). Suppose that the polygon F has 
n vertices. Then p(A, Up) > t„, where 

h = 

h =0.031 

t 5 = 0.076 

f 6 = 0.121 

t 7 = 0.166 

t» = n > 8. 

The proof of this theorem is rather long. The proof extends for twenty pages in 
ifTBI pp. 19-38]. The analogous estimate in the proof of the Kepler conjecture takes a 
full thirty pages ifTSl pp. 126-156]. We cannot pretend to give justice to the proof under 
the page constraints imposed on this version. The reader is referred to the two articles 
just cited for full details of the proof. This article gives a general summary of the ideas 
of the proof, with references for the reader who wishes to pursue the proof in greater 
detail. 

6.1 Verifications in low dimension 

The first two cases n = 3,4 of the theorem can be handled directly with interval arith- 
metic, because they are explicit nonlinear inequalities involving a small number of 
variables. The case n — 3 can be expressed as a nonlinear optimization problem over 
a tetrahedron whose edge lengths vary in length between 2 and 2td d- In other words, 
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it is a minimization problem on the six-dimensional domain [2,2t dod ] 6 . This is readily 
treated by interval arithmetic |26|. 

The case n — 4 can also be directly proved with interval arithmetic. Here the opti- 
mization runs over a nine-dimensional domain. The quadrilateral face F = (vi, V2, V3, V4) 
is parameterized up to rigid motion by the nine variables (3 coordinates for each of four 
points minus the 3 dimensional group of rotations). Monotonicity arguments reduce 
the configuration to a seven-dimensional domain. (Two of the points v, can be rescaled 
v, 1— > Ai'j with < X < 1 until a constraint is met, because parallel shifts in faces of 
a truncated Voronoi cell towards the origin are decreasing in volume.) The inequal- 
ity ju(A, U) > ?4 on a seven-dimensional domain can be proved directly by interval 
arithmetic |26|. 

6.2 Strategy: superadditivity 

Define constants D dod (n, k) for 

n > 3, < k < n, n + k > 4, 

by 

{0.0155 (/a) = (3,1), 

\t n+k - D dod (3, \)k otherwise. 
For (/ii, k\), (112, k 2 ) in this domain, the following superadditivity holds: 

D d0 d(ni,k\) + D dod (n 2 ,k 2 ) > D dod (n x + n 2 - 2,k x + k 2 - 2). (15) 

In fact, this inequality follows immediately from the definitions and the easily verified 
inequality, for m, n > 4, 

t m + t n > f m+ „_4 + 2D dod (3, 1). 



Note that there are only finitely many cases involved in the verification of this identity, 
because t„ is constant for n > 8. 

One of the basic strategies of the proof is to give a partial triangulation of the face 
F (with n sides) into smaller polygons F\,...,F r . The polygon will have «, sides. 
Let kj be the number of edges {u, v} of F, with \u - v\ > 2t dod . Drawing a diagonal 
increases the number of oriented edges by two, and this is the reason for the shift by 



two on the right hand side of (15 1. The proof defines a decomposition of Uf into 
smaller components U-, corresponding to each Fj, gives a bound jt/(A, £/,) > D dod (rii, ki) 
and uses superadditivity <JT3J> to prove the identities: 



s 



r 

fi(A, Ud > V D dod {n h kd > D dod (n, 0) = t n . (17) 



1=1 ;=i 
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The idea is that the objects C7; are lower-dimensional objects than Uf (that is, the 
polygons have fewer edges). The dimension controls the complexity of the estimates. 
Thus a series of inequalities p{A, {/,) > D dod (nj, £,) can be expected to be easier to 
prove than a single inequality p(h, Up) > t„ in higher dimension. 

On the other hand, if the edges of the polygons F, are allowed to get too long, 
numerical experiments show that function /i(A, i/,) tends to become numerically un- 
stable. This prevents an overly aggressive triangulation of F. These experiments lead 
to a restriction of at most 3.2 on the edge lengths. 

A triple (u, v, w) is called unstable if u, v, w are distinct vertices of A* such that 

M < 2t dod , \v\ < 2t dod , \u-v\< 2t dod , |v - w\ < 2t dod , and \u-w\> V8. (18) 

A pair {u, w] is unstable if there exists v such that (u, v, w) is unstable. Otherwise it 
is said to be stable. Unstable edges \u, w) create numerical instabilities and are best 
avoided. 

Remark 6.2.1 (strict inequalities). In most places in the proof, one can be sloppy about 
whether weak or strict inequalities are used, but not here. It is significant that a weak 
inequality \u — v\ > 2t dod is used in the definition of the constant kj. (Deformation 
arguments will be used to decrease \u — v| on a closed interval and the weak inequality 
will keep ki constant.) The strict inequalities are also significant in the definition of 
unstable triples. (Deformation arguments will be used to increase \u\, \v\ on the interval 
[2, 2t dod ] and when the upper endpoint is reached the triple becomes stable.) 



6.3 Construction of subcomponents 

Let F be a face of the graph G(A) and let Uf be the corresponding standard component. 
Represent F as a cycle (vi„... ,v„) with v, € A*. The function p(A, Uf) depends on 
A* only through v\,..,,v n . Thus, for the purpose of the proof of Theorem |6.1| assume 
without loss of generality that A* = {vi, . . . , v„}. 

Say that u e A* is visible from v e A* if {0, u, v} is not a collinear set and if 
aff^ (0, \u, v)) c U F . When this occurs, call the pair {u, v) internal. When {u, v) is inter- 
nal, if the edge \u, v) is added to the graph G(A), the graph continues to be spherical. 

Define {w, v} to be a distinguished pair in U if 

1. \u-v\< 3.2, 

2. {u, v} is internal. 

3. {u, v] is stable. 

Inductively, build a set X of distinguished edges as follows. Start with X = 0. Order 
the distinguished pairs \u, v) by increasing length \u—v\. Considering each distinguished 
edge {u, v] in turn, if it satisfies the non-crossing condition 

aff°(0, [u, v}) n aff°(0, {u , v')) = 0, for {u , v') € X, 

then add it to X, 
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Let G'(A) be the graph on vertex set A* = fvi, . . . , v„} obtained by adding the edges 
X. By the non-crossing conditions, G'(A) is a spherical graph. By the Jordan curve the- 
orem for polygons, Y{A) has two connected components, Uf and the complementary 
region Up>- Each connected component U' c Y(G') is either a subset of Up or Uf. 
Since all the added edges are internal to Uf, there is exactly one component of Y(G') 
that lies in Uf, and that component is equal to Uf. Write [Y(G')Y = [Y(G')] \ {Uf} 
for the set of components of Y(G') internal to Uf. 

Enumerate them U\, . . . , U ,-. Define cj(A, U,) and p(A, Ui), as usual, by dT3}. Then 



( 14 1 and ( 16 1 hold. In fact, the justification given for ( 14 1 holds verbatim in this more 
general context. Let n, be the number of edges of the face F, of G' corresponding to 
Uj. Let k, be the number of edges {u, v) of F, such that \u-v\ > 2t dod . (These are edges 
of G' that do not belong to G(A) and edges of G(A) that have length exactly 2t dod .) 

The function p(A, Uj) depends on A only through the vertices of A on Fj. Thus, 
for purposes of estimating p{A, Ui) for fixed i, assume that A is equal to the set of ver- 
tices of Fj. The estimates can then be expressed locally. This motivates the following 
definition. 

Definition 6.3.1. Let Abe a packing such that € A = A(0, 2t dod ). Let G be a graph 
on vertex set A* consisting of a single cycle containing n > 3 vertices. Suppose that G 
is spherical. Let U € [T(G)] be a connected component ofY(G). The triple (A, G, U) 
is called a local configuration if the following conditions hold: 

1. Every edge {u, v} of G satisfies \u — v\ < 3.2. 

2. If {u, v) is internal in U, then \u — v\ > V& 

3. If{u, v] is internal in U and stable, then \u — v\ > 3.2. 

Let n = n{A, G, U) be the cardinality of A*. Equivalently, n is the number of edges 
in the graph G. Let k = k(A, G, U) < n be the number of edges {u, v} of G such that 
\u-v\> 2t dod . 



The main estimate (Theorem 6. 1 1 now follows from the following refined version 



of the estimate and superadditivity. 

Theorem 6.2. Let (A,G, U) be any local configuration. Let n — n(A,G, U) and k 
k(A, G, U) be the corresponding constants. Assume that n + k > 4. Then 

p(A, U) > D dod (n, k). 

6.4 Deformations 



The proof of Theorem 6.2 is a total induction argument on the cardinality n of A*. 



The induction base case is vacuous, if the induction starts at n = 2, since every local 



configuration has n > 3. Take n > 3 and assume that Theorem 6.2 holds for any local 
configuration with cardinality less than n. 

The strategy of the proof is to deform the local configuration (A, G, U) by moving 
a single vertex v e A* at a time in a way that preserves the constraint of being a 
local configuration, preserves n, is non-increasing in ^(A, U), and is non-decreasing 
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in D c iod(n, k). Under these conditions, any counterexample to the theorem propagates 
to a new counterexample under the deformation. Note that it is easily checked that 
Ddod(n,k) < Ddod(n,k + 1), for all n > 3 and all n > k > 0. Thus, the condition that 
Ddod( n , k) is non-decreasing can be replaced with the constraint that k is non-decreasing 
under deformation. 

The azimuth angle azim(f/, v) is defined for each v e A*. Call v concave in U 
if azim(f/, v) > n. Otherwise, say that v is convex in U. If every vertex v is convex 
in U, then U is a convex set. When U is convex, it is known that U is contained in 
some open half-space whose bounding plane contains the origin. Moreover, if U is 
convex (and n > 3), it is known that u is visible from v in U for any two nonadjacent 
vertices u, v e A*. Note that when U is convex, the conditions on local configurations 
require that \u - v\ > V8, for any two non-adjacent vertices u, v e A*, because {u, v) is 
automatically internal. 



6.5 Deformation at concave vertices 



The most challenging part of Theorem 6.2 is the proof when the local configuration has 



a concave vertex. This subsection sketches the proof in that case. 

The following subsections describe several different deformations. For each, we 
describe the deformation, the starting and halting conditions on the deformation. We 
show that the (A,G, U) remains a local configuration throughout the deformation, that 
/i(A, U) in non-increasing, and k is non-decreasing. 

Let v be concave in U. Let u,w be the two vertices of A* adjacent to v in G. If 
\v-u\ - |v— w\, then the deformation is defined as the continuous motion of v preserving 
|v|, moving along the bisecting plane of {u, w) and increasing |v — u\. If |v - u\ < |v - w\, 
then the deformation is defined as the continuous motion of v preserving |v| and |v - w\ 
and increasing |v - u\. 

The deformation must halt if any of the following conditions are met. (If the initial 
configuration satisfies any of these conditions, no deformation at v occurs.) 

1. For some v' e A* \ \u, v, w], {v, v'} is internal and |v — v'\ < yS. 

2. For some V € A* \ \u, v, w), {v, v'} is internal, stable, and |v - v'| < 3.2. 

3. |v| > 2.2 and \u - v| = |v - w\ = 3.2. 

4. |v| < 2.2, and \u - v\ > 3.07, |v -w\> 3.07. 

By a calculation of derivatives with interval arithmetic, the deformation is non- 
increasing in yu(A, U) [16, Lemma 7.7]. (Although, the function fi(A, U) is potentially 
a function of a large number of variables, the derivative of fi along the deformation 
depends only on the six edge lengths of the simplex {0, u, v, w). This derivative calcu- 
lation is within the reach of interval methods.) The deformation is non-decreasing in k, 
because the length |v - u\ is increasing. 

The deformation preserves |v|, so that the constraint e A = K{Q,2tdod) is pre- 
served. For A to remain a packing, the condition \u - v\ > 2, for u, v e A must hold. 
This article does not repeat the rather technical proof that the condition \u - v| > 2 is 
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preserved, for u, v e A*. The proof runs a couple of pages [16 Lemma 7.6]. It is the 
reparametrization of 1 18, Lemma 12.20]. 

The next constraint is that the deformation should preserve the condition that G is 
spherical. If not, then the deformation produces a situation where v e aff + (0, \u\, u^S) 
for some fixed edge \u\,u-i) of G; or U\ e aff^O, {v, u}) for some fixed U\ of A*. 
The first case is ruled out by IfTBI Remark p. 22], which is the reparametrization of 
lfl8l §12.7,p.l32]. The second case is ruled out by the argument of IfTBI p.27], a 
reparametrization of [18 §12.8,p.l34]. 

The cardinality n of the vertex set A* is preserved. The set of edges of the graph of 
G is combinatorially determined and remains fixed under deformation. In particular, G 
remains a single cycle. Since G is spherical and consists of a single cycle, the set Y(G) 
has two connected components throughout the deformation. The component U evolves 
continuously under deformation. 

The condition for a pair {vi, V2} to be internal is not constant under deformation. 
Nevertheless, IfTBI p.23] (or Ifl8l p.132]) shows that a pair {vi,V2} cannot switch to 
internal when |vi - V2I < V8. When (vi, V2) is stable, it cannot switch to internal when 
I vi - V2\ < 3.2. The enumerated conditions on the internal pairs {vi, V2} in the definition 
of local configurations now follow from the halting conditions. (In fact, the halting 
conditions on internal edges can be replaced with equality, because of the constraints 
on local configurations.) 

The preceding arguments fully justify that the deformation preserves the property 
of being a local configuration, and that any counterexample to the lemma is propagated 
under the deformation. 

There is no loss in generality to assume that the first two halting conditions are 
never met. Indeed, these conditions allow a new stable internal edge {v, v'} to be 
formed. The graph G can be extended to a spherical graph G' by adding the internal 
edge. The component U is partitioned into a disjoint union of two components U\, U2 
of Y(G') and the separating set aff^O, {v, v'}). The vertex set A is the union of A; U A2 
with Ai n A2 = {0, v, V) with corresponding cycles G\ and G2. Both (Ai, G\, U\) and 
(A2, G2, U2) are local configurations. Moreover, yu(A, U) = n(A\, U\) + p(A2, 1/2). By 
the induction hypothesis and superadditivity, the theorem follows in this case. As the 
proof in this case now complete, the following arguments assume that the first two 
halting conditions are never met. 

Thus, the halting condition on v simplifies to 

{\u — v\ — \v — w\ — 3.2, when Ivl > 2.2 

111. 11 
\u - v| > 3.07, |v - w\ > 3.07, otherwise. 

After repeating the deformation at all concave vertices, assume that the halting condi- 
tion holds for each concave vertex. Note that the halting condition at v is incompatible 
with the length conditions in the definition of an unstable triple (v, w\, W2). It follows 
that every internal edge {u, w) at a concave vertex u is stable. By the definition of local 
configuration, this implies that \u — w\ > 3.2. Thus, the hypotheses in the following 
lemma are fulfilled. 

Lemma 6.5.1. Let (A, G, U) be a local configuration with at least one concave vertex. 



Suppose that condition {19 \ holds at each concave vertex. Suppose further \u — v\ > 3.2 
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whenever {u, v} is an internal pair with v e A* concave and u € A* nof adjacent to 
v. Let n,k be the parameters attached to (A,G,U). Finally, assume the induction 



hypothesis that Theorem 6.2 holds for all n' < n. Then 

p(A, U) > D dod (n, k). 

Proof. (Sketch) By definitions, p(A dod ) = f 8 > t n+k - D dod (3, \ )k - D dod (n,k), so it 
is enough to prove p(A, U) > p(A dod ). Define ifr(v, A) to be the angle opposite A in a 
triangle with sides |v|, A, t dod . Recall the right-circular cone rcone(0, v, cos ip(v, A)) from 
the discussion of Tarski arithmetic. The proof breaks into two cases: there are at least 
two concave vertices, and there is exactly one concave vertex. 

Suppose that there are at least two concave vertices. Pick two Vi,V2- Partition 
U into three components U, — U n rcone(0, v,, cos i^(v,, 3.07/2)), for i = 1,2; and 
U = U\(UiU U 2 ). It is known that U\ is disjoint from U 2 [16, Lemma 3.7]. A study 
of the geometry of U\ and U 2 shows the shape of Q. rrmc (A, 0) n t/, and the solid angle 
of Uj depend only on two parameters: |v,| and azim({/, v,). Interval arithmetic gives the 
estimates 

p(A,Ui)>fi(A dod )/2, i = 1,2 

for these two-dimensional objects lfl6l §7.2.6]. On the remaining piece, ju(A, Uo) > 
holds [ 1 8 , p . 1 3 8 ] . The sum of these terms is 



MA !/) = ^/x(A, U i )>p(A dod ). 



i=i 

Now suppose that there is exactly one concave vertex v. In this case, if u e A* \ {v} 
is not adjacent to v, then {u, v} is internal [ 18, p. 140]. This implies that|M-v| > 3.07 for 
all u e A* \ {v}. Consider the deformation that rescales v to decrease its norm v i-> sv, 
for 2/|v| < s < 1. Under this deformation, Q„„ J1C (A, 0) n U decreases in volume and the 
solid angle is unchanged, so that /i(A, U) decreases. Combine this deformation with 
the deformation for concave vertices given at the beginning of this subsection so that 
the constraints in the hypothesis of the lemma are preserved. As before, the induction 



hypothesis is used to avoid the first two conditions of (6.5 1. The constants n,k are 
unchanged; A remains a packing, and so forth. The deformation continue^] until the 
halting condition |v| < 2.2 is satisfied. 

Let Ui = U n rcone(0, v, cos ifr(y, 3.07 - t dod )) and Uq = U \ U\. A study of the 
geometry of Q. trmc (A, 0) n U\ shows that its volume and solid angle only depend on 
two parameters |v| and azim(C/i,v). An interval arithmetic calculation over this two- 
dimensional space, using |v| < 2.2, gives 

H(A, Ui) > p.(A dod ). 

(See H6l §7.2.6]. The constant 1.94159 there is a typo. It should be 3.07 - t dod . The 
typo does not affect the proof.) The inequality /i(A, Uq) > holds for the same reason 



provided in the case of two convex vertices. This completes the proof of Lemma 6.5. 1 



3 A typo in 1 16| incorrectly states |v| = 2. 
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6.6 Deformation at convex vertices 



The results of the previous subsection reduce the proof of Theorem 6.2 to the case 
where U is convex at every vertex. A convex spherical polygon on a unit sphere has 
perimeter at most 2n. A polygon F in G(A) projects to a spherical polygon on the 
unit sphere S 2 . An edge {u,v} of F satisfies bounds |m|,|v| e [2,2t dod ], \u — v\ > 2. 
This implies that every edge of the spherical polygon has arc length at least 8 = 
2 arcsin(l/(2frf„rf)), and that the number of edges is at most seven (2n/6 < 8). 

In this subsection, the geometry is much more explicit than in the previous sub- 
section, because U is convex and F has at most seven sides. Let (G, A, U) be a local 
configuration. Let (n,k) be the associated parameters. This section describes the proof 
of Theorem 6.2 under the total induction hypothesis on n and assuming the truth for 
local configurations already treated. The method, again, is to produce a deformation of 
the local configuration by moving one vertex v at a time. 

Let (A, G, U) be a local configuration with U convex. Let v e A*. Let u, w be the 
two vertices of A* adjacent to v. By convexity, every pair {u, v'}, with v' e A* \ \u, v, w), 
is internal. 



6.6.1 First convex deformation 

Consider the deformation that fixes |v| and |v - w| and moves v to decrease |v - u\. The 
deformation halts (or never starts) once any of the following conditions holds. 

1. For some v' e A* \ {«, v, w), |v - v'\ < V8. 

2. For some V e A* \ \u, v, w), the pair {v, v'} is stable and |v — v'| < 3.2. 

3. azim(t/, v) > n. 

4. There exists V + v such that (u, v', w) is an unstable triple and \u — w\ < 3.2. 

5. |v - w| = VI. 

6. \v-u\ = 2t dod . 

7. \v-u\ = 2. 

As with deformations at nonconvex vertices, the deformation of a local configu- 
ration remains a local configuration. Here, in the convex situation, the proof is more 
elementary, because the geometry is explicit. For instance, the condition that A remains 
a packing follows immediately from the halting conditions, because every nonadjacent 
vertex gives an internal edge. The function fi(A, U) is non-increasing under the defor- 
mation by ifTBI Lemma 7.8]. The value of k in non-decreasing by the halting conditions. 

6.6.2 Second convex deformation 

Let (A, G, U) be a local configuration with U convex. Consider the deformation that 
fixes \u - v| and \u - w\ and moves v to increase or decrease |v|. The direction of the 
deformation is chosen to decrease fi(A, U). The deformation halts (or never starts) once 
any of the following conditions holds. 
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1. For some v' e A* \ {u, v, w}, \v - v'\ < V8. 

2. For some v' e A* \ \u, v, w}, the pair {v, v'} is stable, and |v - v'| < 3.2. 

3. azim(t/, v) > 7r. 

4. |m-v|£{2,2W/, V8). 

5. |w-w| t [2,2t dod , V8). 

6. |v| = 2. 

7. |v| = 2f . 

By an interval arithmetic calculation of derivatives, the function //(A, U) does not 
have a local minimum, provided none of the halting conditions hold iPTBI Lemma 7.10]. 
That is, the deformation can always continue to decrease //(A, U) until a halting condi- 
tion is met. 



6.6.3 Completion of the proof 



Proof. (Sketch) With these two deformations at hand, the proof of Theorem 6.2 can be 
completed. Let n be the cardinality of A*. If n — 3, the inequality of the theorem is an 
inequality in six variables and can be verified directly by interval arithmetic E6l . ifTBI 
§7.4.1]. Now assume that n > 3. Furthermore, by previous estimates, n < 7. 

By induction, it may be assumed that the theorem is established for all n' < n. By 
previous arguments, it may be assumed that the theorem is known for U with a concave 
v (and the same value of n). 

By the induction argument and the reduction to the convex case, there is no loss 
in generality to assume that the first three halting conditions (for both deformations) 
never occur. 

The halting condition (4) of the first convex deformation is rather strange: There 
exists v' + v such that (u, v', w) is an unstable triple and \u — w\ < 3.2. (It was needed in 
the interval arithmetic verifications that prove the monotonicity of //(A, £/)■) Note that 
v and V are both adjacent to u, w when this halting condition holds. This implies that 
n = 4. 

Consider the case n — 4. The dimension of a general configuration is nine, param- 
eterized by four lengths |v| for v e A*, four lengths \u - v] for edges {u, v] of G(A), and 
\u - v] for one internal pair {u, v). Even if the halting condition (4) becomes binding 
at v, deformations can continue at the three other vertices, until some halting condi- 
tion holds at each vertex. Eventually the deformations reduce the dimension of the 
configuration to at most three. Interval arithmetic finishes this case |16, §7.4.2]. 

With the case n — 4 out of the way, the the halting condition for the first convex 
deformation reduces to 

\u~ v| e {2,2t dod , or V8). 

The first convex deformation can be applied at each vertex so that this condition holds 
for every edge {u, v] of G(A). Then, the halting conditions (4) and (5) of the second 
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convex deformation now never occur. The second convex deformation can be applied 
until |v| = 2 or |v] = for each v e A*. 

The local configuration (A, G, U) is now a low-dimensional object. The only re- 
maining continuous parameters are the lengths of (« - 3) internal pairs {u, v) needed to 
triangulate the n-gon G(A). Thus, it has dimension m = n - 3, for 5 < n < 7. 

Unfortunately, the proof does not end here with a simple interval arithmetic calcu- 
lation in low dimension. It does not end here because there is no control on the lengths 
of the triangulating diagonals, and without any such control the calculations are simply 
too numerically unstable. 

A different strategy completes the proof, based on truncated corner cells. Although 
the dimension of this problem is now small, this argument requires several pages. See 
|[T6l pp.30-38]. It is modeled on a published 14-page argument in the solution to the 
sphere packing problem lfT8l §§13.2-13.11], following the same strategy. The final 
parts of this section give a brief summary of the two principal methods that are used. 



6.6.4 Dealing with unstable edges 

If (w, v, w) is an unstable triple, then the halting conditions force \u\ — \w\ — 



\v - w\ = 2. (This relies on the strictness of the inequalities ( 18 i defining stability, as 



mentioned in Remark 6.2.1 ) The value of |v| can be 2 or ltdod- In this final stage of 
the proof, contrary to the constraints of Section |6.3| on distinguished pairs, it is now 
permitted to split the region U into two pieces Uf, Uf separated by aff^O, {u,w}), 
where the pair {u, w} is unstable and \u — w\ < 3.2. One of these pieces is indexed by a 
triangle F. Because of instability, the usual inequality fj.(A, Uf) > Dd ot j(3, 1) does not 
hold. To compensate, stronger inequalities are proved for /i(A, Uf). The deformations 
can continue on the component Uf. Nevertheless, it must be remembered that the 
induction hypothesis and the reduction to the convex case do not cover the stronger 
inequality for ^f(A, Uf) that is now needed. 



6.6.5 Truncated corner cells 

If U has no unstable internal \u, w), then the following argument gives the desired 
bound. The component U is partitioned into n + 1 parts, one U v for each vertex v e A* 
and a final part for the remainder Uq = U\ (U ve v U v ). The function fi(A, U) is a sum of 
terms //(A, U v ) and yu(A, Uo) > 0. The function jj(A, U v ) is a function of the six edges 
of {0, v, m, w] (with m, w adjacent to v) and most of these edges are fixed in length by 
the deformations. So the function fi(A, U v ) is readily bounded with interval arithmetic. 
Each part 

f2 franc (A,0)nc/ v 

is called a truncated corner cell. The defining conditions for U v are 

U v = Un rcone(0, v, cos i/r(v, 1 .6)) Pi H(Q, v, u) n H(Q, v, w), 



where if/(v, A) is the angle defined in Section 6.5 and H(Q, v, v') is the open half-space 
containing v, bounded by the plane through and through the circumcenter of {0, v, v'}, 
orthogonal to the plane of {0, v, v'}. The sets U v , as v ranges over A*, are disjoint from 
one another. We refer the reader to the unabridged version of the proof for details. □ 
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7 Classification of Tame Hypermaps 



This section turns to the problem of classifying a large finite collection of planar 
graphs. For combinatorial simplicity, this classification is phrased in terms of hyper- 
maps, which are defined in the first subsection. The next subsection shows how a 
sphere packing A gives a hypermap. The rest of the section is devoted to the classifica- 
tion problem. A final subsection shows how a counterexample A to the Dodecahedral 
conjecture gives one of the hypermaps classified in this section. 

7.1 Hypermap 

A hypermap is a tuple (D, e, n, f), where D is a finite set, and e, n, / are three permu- 
tations on that set that compose to the identity: e o n o / = /. The elements of D are 
called darts. The permutations e,n,f are called the edge permutation, node permu- 
tation, and face permutation, respectively. (A hypermap was previously defined as a 
finite set D with two permutations /, n, which amounts to the same thing, since e is 
uniquely determined by /, n.) 

If m is any permutation on D, write D/m for the set of orbits in D under m. Simi- 
larly, if G is any group of permutations on D, write D/G for the set of orbits of D under 
G. In particular, D/(e,n,f) is the set of orbits under the group generated by e, n, f. An 
orbit of D under / (n, or e) is called a face (resp. node, or edge). 

A planar graph gives a hypermap by the following procedure. Starting with a planar 
graph, place a dart at each angle (or equivalently at the tail of each directed edge). That 
is, at a vertex of degree k, place k darts, one between each consecutive pair of edges. 
The face permutation has a cycle for each face of the planar graph and traverses the 
darts in a counterclockwise direction around each face. The node permutation has a 
cycle for each vertex and traverses the darts in a counterclockwise direction around 
each vertex. The edge permutation is defined by the relation e o n o / = /. It can be 
interpreted as an involution that pairs a dart associated with the tail of a directed edge 
with a dart associated with the tail of the oppositely directed edge.. See Figure [5] 




Figure 5: A planar graph as hypermap, with faces and nodes 

Hypermaps are the primary combinatorial object used by Gonthier in the formaliza- 
tion of the Four-Color theorem in Coq [9|. Hypermaps, by being purely combinatorial, 
are more convenient to represent on a computer than planar graphs. 
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Not all hypermaps arise from a planar graph in this way. Those that do have two 
special properties. They are involutive and planar in the following sense. The definition 
of planar hypermap is the standard condition on the Euler characteristic, translated into 
the language of hypermaps. 

Definition 7.1.1. 

• The hypermap (D, e, n, f) is involutive, if e is an involution: 

e 2 =I. 

• The hypermap (D, e, n, f) is planar, if 

#(D/e) + #(D/n) + #(D/f) = #D + 2#(D/{e, n, /». 
where #X denotes the cardinality ofX. 

7.2 Packings and hypermaps 

Let A be a packing satisfying e A = A(0,2/^ o j). The graph G(A) is planar. Assume 
that G(A) is biconnected. This subsection describes in greater detail the hypermap 
H(A) attached to G(A). 
For each v e A*, let 

E(v) = \u : \v, u\ is an edge in the graph G(A)}. 

For each u e E{v), there is a half-plane P+(v, u) containing u, bounded by the line 
through {0, v). There is a cyclic order on the half-planes P+(v, u), moving in a counter- 
clockwise circle around the ray emanating from through v. Write cr„ for the cyclic 
permutation on E(v), given by this ordering. 
Define the set of darts by 

D = {(0, v, u, o-yii) : u e E(v)}. 

Define face, edge, and node permutations on D by 

f(0,v,w,u) = (0, w,cr~ l v, v), 
e(0, v, w, u) = (0, w, v, cr w v), 
n(0, v, w, u) = (0, v, u. cr v u). 

A formal calculation shows that (D, e, n, f) is an involutive hypermap. 

The nodes of (D, e, n, f) are in bijection with A* under the correspondence: 

(0, v,w,«)hv£ A*. 

The edges of (D,e,n,f) are in bijection with the edges of G(A) under the correspon- 
dence: 

(0,v,w,u) i-> {v, w). 
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The faces of (D,e,n,f) are in bijection with the faces of G(A) under the correspon- 
dence: 

x = (0, v, w, u) i-» F = ((f-'xh, X3, (fx) 3 , (fxh, . . .), 

where yj, is the third component of the four-tuple y and each face of G(A) is represented 
as usual as cycle (vi, . . ., v„). The set of darts is in bijection with the set of oriented 
edges of G(A) under the correspondence: 

(0, v, w, u) i— > (v, w). 

The graph G(A) is connected. This implies that (e, n, /) acts transitively on D. 

The graph G(A) is planar. If V,E,F are the number of vertices, edges, and faces, 
then V - E + F — 2; or equivalently, V + E + F- 2 + 2E. Under the bijections just 
described, this implies that the cardinalities of these sets satisfy 

#(D/e) + #(D/n) + #(£>//) = #D + 2#(D/(e, n, /». 

Thus, the hypermap (D, e, n, f) is planar. 

Remark 7.2.1. In 4761/ . the basic combinatorial structure is called a planar map rather 
than hypermap. In that article, the combinatorial structure is represented in computer 
code as a finite set of faces 

{F 1 ,F 2 ,...,F r }, 

and each face is represented as a cycle (vi,...,v„) of vertices. This is essentially 
equivalent to a hypermap. This representation is converted to a hypermap by sending 
(v\ , V2, . . . , v„) to the dart (0, vi , vi, o~ Vi Vj). There are n-choices of which vertex v,- to 
list first in the cycle, and by taking all choices, n darts are obtained. Running through 
all faces in this way, all darts are constructed. In the opposite direction, an earlier 
argument describes how a face of the hypermap gives a face ofG(A), expressed as a 
cycle. 

7.3 Tameness 

This subsection defines a collection of hypermaps called tame Voronoi hypermaps. 
The classification of these hypermaps, up to isomorphism, is one of the main steps of 
the proof of the Dodecahedral conjecture. Bauer and Nipkow have formally proved 
the classification in Isabelle ll29l . This Isabelle proof was originally designed for the 
classification of the sphere packing problem. Nipkow and McLaughlin have modified 
that proof to cover the Dodecahedral conjecture as well. The modified Isabelle proof is 
found at ll26l . 

Let H = (D, e, n, /) be a hypermap. A face of H (that is, an orbit of D under f) is 
said to be a triangle, quadrilateral, pentagon, etc. if the cardinality of the orbit is 3, 4, 
5, respectively. Two nodes are said to be adjacent if there is an edge \x,y} of H such 
that x belongs to one of the nodes and y belongs to the other. 

Let v be a node of H (that is, an orbit of D under «). A node v is said to have 
type (p, q, r), if the cardinality of v is p + q + r and if there are p triangular faces, q 
quadrilateral faces, and r other faces that share a dart with v. The cardinality of a node 
is also called its degree. 
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* 


* 



Table 1 : Vertex types 



Define constants b(p, q) by Table [T] If (p, q) falls outside this table, or if the entry 
is marked *, then set b(p,q) = n(Ad d)- Let t„, n > 3, be the collection of constants 
defined in Theorem l6.ll 

A weight assignment of a hypermap (D, e, n, /) is a function w : D — > R that is 
constant on faces: w(fx) = w(x) for x e D. A weight assignment is said to be a 
Voronoi weight assignment if the following properties hold: 

1 . If the face containing x has cardinality to, then w(x) > t m . (In particular, w(x) > 
for all jc.) 

2. Let F c D be any face with cardinality m > 5. Let y € F. Let V be a set of 
nodes, each meeting F, such that no two are adjacent to one another. Assume 
that the type of each node of V is (4, 0, 1). Let X = (IIV) \ F; that is, the set of 
darts in nodes in V except those in F. Let m' be the cardinality of V. Then 

w(y) + V w(x) > t„, + 0.016 m'. 

xeX 

3. If the node of x has type (p, q, 0) and degree m — p + q, then 

m 

^ w(n'x) > b(p,q). 
i=i 

The total weight of a weight assignment w is defined to be 

■re [£>//] 

where [£>//] is a set of representatives of the orbits of D under /. 

Definition 7.3.1. H — (D, e, n, f) is said to be a tame Voronoi hypermap if the following 
conditions hold. 

1. H is an involutive, planar hypermap. 
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2. H is connected. That is, D is a single orbit under (e, n, f). 

3. (Simple face) Every face of H meets every node of H in at most dart. 

4. The number of nodes is at least 13. 

5. The cardinality of each face of H is at least 3 and at most 7. 

6. (Triangle types) Let Vi, V2, V3 be any three nodes of H such that v,- is adjacent to 
Vjfor each i + j. Then there is a dart x,- € V; such that {xi, X2, x-j] is a face of H. 

7. There are never two nodes of type (4, 0, 0) that are adjacent to one another. 

8. ( Quadrilateral types ) Let V\ , V2 , V3 , V4 be any four distinct nodes of H such that 
Vi is adjacent to V;+i for i — 1,2,3,4 (setting V5 = v\). Then darts X; e V; can be 
chosen so that the faces containing Xifall into one of the four patterns depicted 
in Figure^ 

9. The cardinality (degree) of a node of type (p, q, r) is at most five if r > 0. 

10. The degree of each node of H is at least 2 and at most 6. 

11. There exists a Voronoi weight assignment of total weight at most ^(A^). 




Figure 6: Quad types 



7.4 Classification 

Two hypermaps (D, e, n, f) and (£)', e', n',f) are properly isomorphic if there is a re- 
jection between D and D' that is equivariant for the face, node, and edge permutations. 
Each hypermap (D, e, n, f) has a mirror image: 

{DJn,n l J- 1 ). 

An improper isomorphism between two hypermaps is a proper isomorphism between 
one hypermap and the mirror image of the other. Two hypermaps are isomorphic if 
there is a proper or improper isomorphism between the two hypermaps. If a hypermap 
is involutive or planar, then so is every isomorphic hypermap. 

There is an archive of tame Voronoi hypermaps |26l . This archive contains 206 
hypermaps. 
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Theorem 7.1. If H is a tame Voronoi hypermap, then it is isomorphic to some hyper- 
map in the archive. 

The proof of this theorem relies on the piece of computer code described in Sec- 
tion 3.4 The reader is referred to that section for a description of the details of the 



theorem. 

7.5 Counterexamples are tame Voronoi hypermaps 

The following theorem proves that every potential counterexample to the Dodecahedral 
conjecture gives a tame Voronoi hypermap. In particular, the classification of tame 
Voronoi hypermaps gives an explicit case enumeration of the possible combinatorial 
structures of a counterexample. The following section will eliminate each case in the 
enumeration. This will eliminate all possible counterexamples to the Dodecahedral 
conjecture. 

Theorem 7.2. Assume that A is a counterexample to the Dodecahedral conjecture. 
Without loss of generality, assume that € A = A(0, ltdod); that the cardinality of 
A* is at least 13; and that G(A) is biconnected. Let H — (D,e,n,f) be the hypermap 
attached to the graph G(A). For every dart x e D in face F, set 

w{x) = p(A, U F ). 

Then H is a tame Voronoi hypermap and w is a Voronoi weight assignment on H of 
total weight at most p(Ad d). 

Since the definition of tame Voronoi hypermap is a long enumeration of different 
properties, the proof of this theorem breaks into a long enumeration of lemmas, each 
establishing one property. The statement of the theorem specifies the weight assign- 
ment w. The verification that w is a Voronoi weight, breaks into separate lemmas for 
each property in the definition. The article [16| devotes many pages to the proofs of 
these lemmas. This articles sketches the proofs and refers the reader to the fuller ver- 
sion for details. Turn to an item by item discussion of the properties. The first several 
are elementary. 

7.5.1 Involutive 



The hypermap H is involutive and planar. This has been established in Section 7.2 



7.5.2 Connected 



The set of darts D is a single orbit under (e, n, /). This is also contained in Section 7.2 
It follows directly from the connectedness of G(A). 

7.5.3 Simple 

Every face of H meets every node of H in at most dart. This is a direct conse- 
quence of the biconnectedness of G(A). Indeed, if by following the face permutation 
x, fx, f 2 x, . . . , fx = y, the darts x + y lie at the same node v; then v is an articulation 
vertex of G(A), and the graph is not biconnected. 
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7.5.4 Node cardinality 



The number of nodes is at least 13. This is a consequence of the assumption that the 
cardinality of A* is at least 13 and the bijection in Section [772] between nodes of the 
hypermap and A* . 



7.5.5 Face cardinality 

The cardinality of each face of H is at least 3 and at most 7. By definition, the face 
map on a dart x takes the form 

x = (0, v, w, u), fx = (0, w, cr~ v, v), f~x = (0, u, v, o- u v). (20) 

In a biconnected graph with more than two vertices, every vertex has degree at least 
two. Thus every node of the hypermap has degree at least two. Thus, o~ u v + v and 
fx + f~ x X. Also, w,u e E(v), which does not contain v. Thus, the form of x, fx, and 
f x in (20 1 shows these darts are distinct, and the face contains at least three distinct 
darts. 



If some face of H has cardinality at least 8, then by Theorem 6. 1 
//(A) = ^ fi(A, U)>t s = yu(A dorf ). 



Ue[Y(A)] 

Thus, A is not a counterexample, as was assumed. 



7.5.6 Triangle types 



Let Vi, \>2, V3 be any three nodes of H such that v,- is adjacent to vjfor each i + j. Then 
there is a dart Xi e v; such that {X\,X2,x^\ is a face of H. This is a restatement of 
Lemma 5.3.2 in terms of the combinatorial properties of hypermaps. 



7.5.7 Adjacent degrees 

There are never two nodes of type (4, 0, 0) that are adjacent to one another. If there are 
two adjacent nodes of type (4, 0, 0), then the graph takes the shape of Figure [7] This is 
an impossible configuration in a packing A for purely geometric reasons. It has nothing 
to do with the value of /u(A) and volumes of truncated Voronoi cells. The impossibility 
proof appears as [16, Lemma 3.8]. It is a reparametrization of |[T3l Prop.4.2]. This is 
one of the most delicate reparametrizations. 



7.5.8 Quadrilateral types 

Let vi, V2, V3, V4 be any four distinct nodes of H such that v,- is adjacent to v,+i for i — 
1,2,3,4 (setting vs — v\). Then darts Xj e v,- can be chosen so that the faces containing 
Xjfall into one of the four patterns depicted in Figure^ Let G' be the spherical graph 
on the vertex set {vi, v%, V3,V4} with edges {v;,V; + i} for i - 1,2,3,4. Then by Jordan 
curve theorem for polygons, Y(G') consists of two connected components. The result 
of ifTBI Lemma 3.8] cited in the previous proof states more precisely that exactly one 
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Figure 7: Two adjacent vertices of type (4, 0, 0) 



connected component U of Y(G') has solid angle less than 2n, and that t/ n A contains 
at most one point. If U n A is empty, then the first pattern of Figure[6]occurs. 

In remaining cases, (/ flA contains a single point vo e A. Again, by the Jordan 
curve theorem and the planarity of G(A), all edges {vo, u) in G(A) have the form {vo, v,-} 
for i = 1,2,3,4. Section 7.5.5 shows that each node has degree at least 2. Thus, vo is 
adjacent to 2, 3, or 4 of the vertices v,-. Figure [6] gives all such connection patterns of 
vq with Vj, except the one shown in Figure [8] 



Figure 8: A quad type that is excluded 

Thus, it is enough to show that pattern of Figure[8]does not occur in any counterex- 
ample to the Dodecahedral conjecture. This graph contains a triangle F and a pentagon 
F' , with corresponding standard components U - Uf and U' = Uf s [T(A)]. An 
important estimate |[T6l Lemma 10.1] gives that 

//(A, U)+n(A, V) > 0.168. 

If there is some other face with n > 4 sides, then Theorem |6. l| gives 

M(A) > fi(A, U)+fi(A, U') + t n > 0.168 + 0.031 > ft(A dod ). 

Otherwise, pick any four vertices in A* other than Vo, Vi, . . . , V4. Let U\,...,U r be 
the standard components indexed by triangles of G(A) that contain one of these four 
vertices. Another estimate |[T6l Lemma 5.2] gives 

r 

V [£(A, Ui) > 4(0.004) and hence fi(A) > 0.168 + 4(0.004) > /u(A dod ). 
This shows that A is not a counterexample. 



41 



7.5.9 Degree 



The degree of a node of type (p, q, r) is at most five if r > 0. Indeed, it is a direct 
consequence of the inequalities of |26| that if U\, . . . , U n , for n — p + q + r > 6, are the 
standard components of the faces around the node, then 

n 

2^//(A, [/,) > fx(A dod ). 

i=l 



7.5.10 Degrees 



The degree of each node of H is at least 2 and at most 6. Section 7.5.5 has already 
shown that the degrees are at least 2. Let (p, q, r) be the type of node v. If r > 0, then 
the previous property bounds the degree at five. Assume r — 0. The proof in this case 
is deferred until Section |T.5.12| 



7.5.11 Total weight 

There exists a Voronoi weight assignment of total weight at most n{Kd Q d). By definition, 
a counterexample A is a packing such that 

/i(A)= Yj KKU) < ni^dod). (21) 

Ue[Y(A)] 



The set [T(A)] of standard components is in bijection with the faces of H. By the 
definition of w in the statement of the theorem, ( pT} can be rewritten as 

^ W(X) < p(Klod)- 

xe[D/f] 

This is exactly what it means for w to have total weight at most p{Adod)- 



7.5.12 Voronoi weight 

The weight assignment w is a Voronoi weight. This can be expressed as the following 
claims. 

1. If the face F has cardinality n > 3 then p(A, Uf) > t„. 

2. Let F c D be any face with cardinality m > 5. Let y € F. Let V be a set of 
vertices of F, such that no two are adjacent to one another. Assume that the type 
of each node of V is (4, 0, 1). Let X be the set of triangles at the vertices V. that 
is, the set of darts in nodes in V except those in F, Let m' be the cardinality of 
V. Then 

//(A, U F ) + V //(A, U F ,) > t m + 0.016 m . (22) 

F'eX 
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3. If the node of x has type (p, q, 0) and degree m = p + q, then 

in 

Y J KKU i )>b(p,q). (23) 



The first of these claims is a direct consequence of Theorem 6.1 Consider the 
second claim. Form the graph G'(A) and partition Uf into subcomponents by the 
algorithm described in Section [63] A result in Tarski arithmetic states that two internal 
pairs \u, w) and \u', w'} with \u— w\, \u'— w'\ < V8 do not cross each other in the sense of 
Section 6.3 [ 16, Sec. 7.2.1]. It follows that every internal pair {u, w) such that \u - w\ < 



V8 forms an edge of the graph G' . The bound on t m is obtained by superadditivity, 
from a collection of inequalities ju(A, {/,-) > Ddod{n,k) for each subcomponent. 

Let v be a node in V. Let (0, v, u, w) be the dart in F at node v. Let U v be the sub- 
component indexed by the triangle {v, u, w) in G' . At a node of type (4, 0, 1), by angle 
considerations, the pair \u, w) is necessarily internal. Let U' v . . . , U' 4 be the standard 
components indexed by the four triangles at v. By summing over V, the following two 
inequalities imply (|22]i. 



4 

J^H(A,U'i) > 0.016. 



1. \f\u-w\ > V8, then 



2. If \u - w\ < V8, then 

4 

MA, U v ) + ^ ^(A, U\) > D dod O, 1) + 0.016. 



;=1 



These inequalities are obtained by summing over interval arithmetic inequalities for 
each term fi(A, •)■ A detailed proof appears at [16 Theorem 8.1]. 

Consider the final claim. The argument at ifTBI p. 14] goes as follows. Let v e A* be 
a vertex. To have type (p,q,Q) means that there are p triangles and q quadrilaterals in 
the graph G(A) at the vertex v. Interval arithmetic can be used to compute lower and 
upper bounds on the azimuth angles azim(f//r, v) when Uf is a triangle or quadrilateral. 
These bounds are IfTBI F.2.1,F.4] 

triangle: 0.856147 <azim(£/,v) < 1.88673, 

quadrilateral: 1.15242 <azim(£/,v) < 3.25887. 

Ill 

^ azim(t/ ; , v) = In. (24) 

By the bounds on the azimuth angles, this equality can only be satisfied for special 
(p, q). Anything else is a geometric impossibility. The feasible pairs {p, q) are listed in 
IfTBI Lemma 6.1]. According to this list, p < 7 and q < 5. 

Interval arithmetic methods establish a list of nonlinear inequalities relating azim(f/f, v) 
to fi(A, Uf) when F is a triangle or quadrilateral [26]. If free variables azim(F) and 
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mu(F) are substituted into these inequalities for azim([/p , v) and p(A, Uf), then the re 
suiting inequalities are linear in these variables. A system of linear inequalities results 
(This is the method of linear relaxation.) A lower bound on the left-hand side of ((23] 
is the solution of the linear program 



mm 

' F 



i^]muCF) 
subject to 

^ azim(F) = In, 

F 

and to the system of linear inequalities. The linear program is run for each p, q and a 
constant b(p, q) slightly smaller than the minimization was picked. This gives the table 
of values. 

If A is a counterexample to the Dodecahedral conjecture with a vertex of type 
(p, q, 0), then 

piAdod) > y"(A) > b(p, q). 

An inspection of the list of constants b(p,q) in the table shows that that this implies 
that p + q < 6. That is, the degree of every vertex of type (p, q, 0) is at most 6. This is 



the property needed in Section 7.5.10 



8 Linear Programs 



This section discusses Theorem 8.1 It is one of the main steps in the proof of the Do- 
decahedral conjecture. The discussion begins with the terminology used in the state- 
ment of the theorem. 

Definition 8.0.1. A hypermap system is a pair (H, <E>), where H = (D,e,n,f) is a 
hypermap, and <D is a finite set constraints on H. More precisely, let V be the vector 
space of real-valued functions on D. Each <p e Q> is a boolean valued function <p : 
V e — > {true, false}/or some (. (Assume ( is independent of $ e <&). 

The hypermap system (H, 4>) is said to be feasible, if there is some x = (xi, . . . , X[) e 
V e such that (p(x) holds for all <p e <£. Otherwise, the system is infeasible. 

In our hypermap systems, the constraints <D are generated from a list of about one 
hundred generic parametrized constraints, with a parameter running over the darts in a 
hypermap. For example, a parameterized constraint c(H, x) might be interpreted as an 
inequality relating fi(A, U f), the solid angle of U f, and azimuth angles, whenever F is 
a triangular face containing the dart x. As x runs over the darts in a hypermap H, the 
generic constraint c(H, x) evaluates to a constraint for each dart x in every triangular 
face of H, yielding dozens of particular constraints for H. We do not list the generic 
constraints here, but they appear in 1 16 1 and also at ll26ll . 

For each tame Voronoi hypermap H, the generic constraints specialize to a par- 
ticular set of constraints <$>dod{H). Each set $>dod{H) contains about 5000 constraints. 
The hypermap systems (H, ^dod(H)) arising in this way are called Voronoi hypermap 
systems. 
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Theorem 8.1. Let H be any tame Voronoi hypennap. Then (H, <t>d d(H)) is infeasible. 



This proof is carried out by computer as a collection of linear programs. This is 
one of the three major parts of the proof of the Dodecahedral conjecture that have been 
carried out by computer. This article describes the relationship between the feasibility 
of (H, <$>) and a linear programming feasibility problem. It also describes some details 
of the implementation of the code. 



Section 7.4 enumerates of all tame Voronoi hypermaps. Thus, the proof of Theo- 
rem 8.1 may proceed case by case. Let's focus attention for a moment on one tame 
Voronoi hypermap H = (D,e,n,f) and the corresponding Voronoi hypermap system 
(H, O). A simple strategy will show that it is infeasible. For some t e N, each con- 
straint <p is a function on V c , where V is the vector space of real-valued functions on D. 
Thus, V c can be identified with R m , where m - (#(D). An inspection of the form of 
the generic constraints <p e O described above reveals they all have a very special form. 
They are all linear constraints on R m . 

Some of the linear constraints carry guard conditions. That is, some constraints 
have the form 

(Ax </?)=> (A'x < b'\ (25) 

for x e R" ! , and various matrices A, A' and vectors b, b'. (The vector inequality a < b 
means that a, < bj for every component of the vectors a,b.) The constraint (Ax < b) 
is called a guard condition. Variations are allowed in which some of the inequalities in 
the guard condition are weak and some of the inequalities in the consequent are strict. 

The collection of all inequalities that do not have a guard condition is a system of 
linear inequalities. Standard linear programming packages can be used to determine 
whether this system of linear inequalities has a feasible solution. If this linear program 
is infeasible, then the hypermap system (H, <t>) is clearly also infeasible. When this 
happens, a proof of the infeasibility of (H, results. 

When this fails, the constraints with guard conditions are used. The introduction of 
a constraint that has a nontrivial guard condition involves multiple steps. The constraint 



( 25 i can be rewritten in logically equivalent form as 



(A { x > bi) V ■ ■ • V (A r x > b r ) V (A'x < b'), 



where A, and bj are the rows of A and b. Taking each disjunct in turn, one linear 
inequality at a time is added to the system of linear inequalities, and the resulting 
system is shown to be infeasible. When each systems are infeasible, then (H, O) itself 
is infeasible. 

This discussion may give the impression that a great many linear programming 
feasibility problems are created in this manner. In practice, nearly all of the hyper- 
map systems are eliminated in the first pass, without requiring recourse to the guard 
conditions. 



8.1 Counterexample implies feasibility 

Theorem 8.2. Let A be a counterexample to the Dodecahedral conjecture. Let (H, <I>) 
be the Voronoi hypermap system attached to A. Then (H, O) is feasible. 
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Proof. It is enough to give some assignment of x = (x\,...,x£) e V such that <p{x) 
holds for all (f> e <t>. Each x,- is a real-valued function on darts. The notation for the 
functions x,- has been set up in a way that suggests the assignment. Let a = (0, v, u, w) 
be a dart. Let F be the face containing a. Make the following settings. The first row of 
this table gives the functions x, in the notation of [26|. The second row gives the real 
number x,(a). 

yn ye sol azim mu omega 

|v| \v-u\ sol(U F ) azim(t/ F ,v) h{K,Uf) w(A, Up ). 

All of the predicates O can be shown to hold for this assignment, either as a con- 
sequence of definitions, as consequences of geometrical facts, or as consequences of 
interval arithmetic calculations. 

The interval arithmetic calculations for the predicates with guard conditions were 
verified in two stages. The polygonal face is triangulated, and the standard component 
Uf is partitioned in a corresponding way into parts U\, . . . , Uk- In the first stage, the 
guard hypotheses were used to obtain interval arithmetic bounds on the lengths of the 
internal edges of the triangulation. In the second stage, these edge length bounds are 
used as hypotheses in further interval arithmetic bounds that give lower bounds on 
<j(A, Ui) > c,. Then 

k k 

w(A, U F ) = J] w < A > U i) > Yi Ci = °- 

i=i i=i 

This becomes a general predicate on darts of the form: 



oraega(o') > c. 



Since no tame Voronoi hypermap is feasible, by the argument of Section 2.2 the 
Dodecahedral conjecture is established. 
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